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Method for non-invasively measuring respiration volume whereby weighting factors are non-invasively determined 
by (a) totaling delta values over a baseline period of substantially steady state breathing for obtaining a rib cage signal, (b) 
totaling delta values over a baseline period of substantially steady state breathing for obtaining an abdominal signal, (c) di- 
viding the average variability of the mean of the total of the delta values for one of either the rib cage or abdomen signals 
by the average variability of the mean of the total of the delta values for the other of either the rib cage or abdomen sig- 
nals; and (d) multiplying the other signal by the quotient derived from step (c). 
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SINGLE POSITION NON-INVASIVE CALIBRATION TECHNIQUE FOR 
RESPIRATION MONITORING APPARATUS 

Technical Field 

TItia invention pectains to methods and apparatus for 
measuring respiration volume and, more particularly, to 
such •♦thods and apparatus which measure respiration vol- 
5 ume by separately measuring and then summing the contribu- 
tions from a plurality of torso portions, such as the rib 
cage ana abdomen. Most particularly, this invention 
10 pertains to a calibration technique for weighting signals 
indicative of the contrmutions from the torso portions 
whereDy the sum of the signals is proportional tc respira- 
tion volume. 
Background Art 

Unitea states Patent No. 4,308,872 of January 5, 
1982, entitled Methoa And Apparatus For Monitoring Respi- 
ration, the contents of which are incorporated herein by 
reference in their entirety, discloses a method and appa- 
ratus for quantitatively measuring respiration volume. 
25; The method disclosed in that patent comprises looping 

first and second extensible conduct: :s scout the rib cage 
and abdomen, separately and simultaneously measuring the 
30 inductances of the conductors dun.-: respiration, 

weighting the measured inductances to reelect the differ- 
ent contributions of tr.e rio cage and accomen to respira- 
tion volume, and summing the weighted measured inductances 
to obtain actual respiration volume. 

as notea, practice, of tne cecnr.i^ue disclosed in tne 
patent requires weighting or calibrating tne mauctances 
40 measured by the aDdomen ana rib cage conauctors. To 
effect calibration it is necessary to determine the 
weighting . factors K and L to satisfy tne following equa- 
45. tion: 

V =• K-RC + L * AB [.EQUATION A) 

where V is total respiration volume, RC is the rib cage 

50 



35 
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1 

I contribution to respiration volume as measured at the rib 
C49* conductor and AB i3 the abdominal contribution as 
measured at the abaominal conductor. Patent No- 4,308,572 
. 3 discloses a specific tecnnique for determining tne values 
for the weighting factors K ana L. 

In accordance with that disclosure, a spirometer is 

employed during the calibration procedure. With the 
10 „ 

patient in a first position, such as standing, a simulta- 
neous set of readings are recorded from the outputs of the 
spirometer, cne no cage conductor, and the abdominal con- 
15 due tor. This ls repeated with the patient- in a second 
position such as supine. At this point, there are two 
sets of vai-es for V, RC, and AB which satisfy Equation A. 
Thus, two equations having two unknowns, the constants K 
and L f may ce -riccen. From these, the weighting factors 
K and L may esteemmed by employing well known 
techniques :o: solving simultaneous equations. Thus: 
K * AB 1 . V 2 - AB 2 . V x 

(EQUATION BJ 
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RC 2 ABj/- RC 1 AB 2 
L = RC X . V 2 - V x . RC 2 



i EQUATION CJ 



AB 2 RC X - AB X RC 2 

35 Tfcs denominators oc Equations o ir.c w rr.ay, ascending 

upon the recorded values, aepcoaen :: -^qjal zero. 
Clearly, wnen this, happens, —e va.-es retained cor K ana 
L will be inaccurate, thereoy sowing any measurement 
based on sucr. weighting factors. TT.js, eacn. time the 
denominators 3 £ Equations 3 and C 2?cr cacn or equal zero, 
a new set of readings must be taxen, -nereoy increasing 
the time required cor calibration. 

In U.S. Patent No. 4, 37 3, 534 of cecr.uary 15, 1983, 
entitlea Mecr.oc And Apparatus For Caiiorating Respiration 
50 Monitoring System, tne contents ot which are also incor- 
porated herein oy reference in their entirety, an alter- 
nate method ana apparatus for a. graphmg-oased technique 
55 for determining tne weighting factors K and L is 

disclosed. As in the simultaneous equation technique of 
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U.S. Patent No. 4,308,872, a spirometer or other device 
for independently measuring respiration volume is employed 
during the calibration procedure. With the subject in a 
first position, readings from the spirometer, the rib cage 
conductor and the abdominal conductor are simultaneously 
recorded for a plurality of breaths, preferably at least 
three in numoer. This is repeated with the subject in a 
second position. For each breath, the rib cage and abdom- 
inal readings are divided by the spirometer reading. That 
is, the values RC/V and AB/V are obtained for each breath, 
15 where V is the respiration volume. as measured by the 

spirometer, RC is the rib cage reading from the uncali- 
brated rib cage conductor, and AB is the abdominal reading 
20 from the uncalibrated abdominal conductor. The points 

(RC/D, AB/V) for each breath are next plotted on a graph 
whose axes are RC/V and AB/V, and a line approximation is 
drawn through tnese points. The line may be drawn by vis- 
ual approximation, although preferaoly it is determined by 
the least squares technique. The .:ne is then extended 
through the x ana y axes. The reciprocals of the x and y 
intercepts oefme tne weighting factors K and L; i.e. the 
reciprocal of tne intercept of tne RC/V axis defines the 
weighting tactor K for the no cage ana the reciprocal of 
35 . the intercept ot the AB/V axis defines tne weighting fac- 
tor L for the abdomen. Pcetsraoiy, all or the foregoing 
calculations are carried out by a rr.icrccc Jcessor or other 
data processor which performs tne calculations and yields 
values for tne weighting factors K una L . 

A drawback of the methods and apparatus disclosed in 
U.S. Patents Nos. 4,308,872 and 4,373,534 is the require- 
ment that sets of data points or values from the abdom- 
inal and rio cage conductors and from the spirometer or 
other respiration volume measurement device be obtained 
50 with the subject for two different distributions of venti- 
lation — i.e. in two separate positions. Where a sub- 
ject's physiological .condition prevents or dictates 
55 against movement from one position to another, however, 
calibration of the rib cage and abdominal conductor 
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^ contributions cannot be readily carried out in accordance 

with- the known methods. 

^Another drawback of these prior art methods and appara- 
* Ac- 
tual is^ that in certain applications, such as neo-natal 

monitoring/ it is not practical to calibrate the apparatus 

using an independent respiration volume measuring device 

such as a spirometer. For example, since the above 

10 calibration techniques require airway connection to a 
spirometer or other similar device, significant time is 
required to carry out trie procedure. This is often 

15 . unacceptable to new-born nursery staff where time is at a 
premium. 

Yet another problem, with these calibration 
techniques and apparatus is that they rely on the 
assumption that all air movement in the respiratory system 
is between the rib cage and spirometer or the abdomen and 
the spirometer. In fact, there also exists movement of 
25 air between the rib cage and the abdomen during normal 
respiration. This RC-AB exchange of air is pendelluft 
that occurs continuously and with varying degree in respi- 
30 ration. The known methods do not incorporate an allowance 
for pendelluft. 



35 



40 



20 



45 



50 



3D 



60 



WO 87/04332 



PCT/US87/00217 



-5- 

5 Another known quantitative calibration technique is described in 
Single Position Tidal Breathing Calibration of the Respiratory 

^ Inductive PI ethy sinograph, Watson et al, Amer. Rev of Respiratory 
Diseases, Vol. 129, p. A256 (1984). According to that 
technique, with the subject in a single position or posture, 

15 readings from a spirometer (SP) , the rib cage conductor (RC) , 
and the abdominal conductor (AB) are simultaneously recorded for 
a predetermined period of preferably at least one full breath. 
The curves (SP, RC) and (SP, AB) are then plotted from the 
recorded data and each resulting curve is closed by a straight 

25 line connecting its beginning and end point. The resulting loop 
areas are then calculated as by integration. Then, using any 
selected data points simultaneously recorded from the spirometer 

30 

(SP) and from the rib cage (RC) and abdominal (AB) conductors, 
weighting factors (for the rib cage scaling amplifier and/or for 
35 the abdomen scaling amplifier) may be determined. This 

calibration technique also suffers from the drawback that it is 
invasive in the sense that it requires a spirometer or other 

40 

device measuring respiration volume at the mouth. 



45 



50 



55 



60 



WO 87/04332 PCT/US87/00217 



-6- 

5 grje£ Desc ription Of The Drawings 

In the drawings, wherein like reference numberals denote 
similar elements throughout the several views: 

10 

FIG. 1 is a diagrammatic representation of a portion of a 
system for non-invasively monitoring respiration volume; 
15 FIG. 2 is a block diagram of a complete system for 

non-invasively monitoring respiration; and 

FIG. 3 is a graphic representation illustrating the delta 

20 

values for the rib cage and abdominal signals. 
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5 Bast Mode For Carrying Out T he Invention 

Referring now to the drawings, apparatus for measuring 
respiration volume of the type disclosed in U.S. Patent No, 

10 

4,306,872 is shown in FIGS* 1 and 2, Two extensible conductive 
loops 12, 14 are secured in any suitable fashion to elastic 
15 tubes 16, 18, respectively, such that the conductors 12 and 14 
extend respectvely about the rib cage and abdomen of the subject 
10. As the subject 10 breathes, the elastic tubes 16, 18 and 

20 

conductive loops 12, 14 expand and contract, resulting in 
changes in the inductances of the loops. After the inductance 
25 of each loop is converted to a proportional signal, the signals 
are calibrated and then summed to provide a signal indicative of 
tidal volume. Calibration of the signals from the rib cage and 

30 

abdomen conductors 12 and 14, respectively, is necessary because 
the relative contributions of the rib cage and abdomen to tidal 
35 volume vary from subject to subject and even in a single subject 
with different postures, e.g. standing, supine, etc. 

Suitable apparatus for converting the inductances of the 

40 

conductors 12 and 14 to proportional electrical signals, 
calibrating those signals to reflect the proper relative 
45 contributions, and then summing those signals to provide a 

signal indicative of tidal volume are known to those of ordinary 
skill in the art. One such apparatus is disclosed in U.S. 

50 

patent no. 4,306,872. Another suitable apparatus is marketed by 
Nims, Inc., Miami Beach, Florida under the model designation 
55 Respigraph ™. Such apparatus are generically illustrated in 
block diagram form in FIG. 2, where the blocks 20, 22 represent, 
respectively, appropriate circuitry for converting the 

60 

inductances of the rib cage and abdominal conductors 12, 14 to 
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5 signal* suitable for further processing. The scaling amplifiers 
24, 26 represent circuitry for calibrating the signals from the 
rib cage and abdomen, respectively, to reflect the relative 

10 

contributions of the rib cage and abdomen to tidal volume. Once 
the scaling amplifiers 24, 26 are properly calibrated, the 
15 resulting signals may be summed, as by the summing- amplifier 23 
in FIG. 2, to yield a signal indicative of tidal volume. The 
output signal from the summing amplifier 28 as well as the 

20 

output signals from the two scaling amplifiers 24 and 26 may 
then be displayed as on a graphic recorder 30 or a digital 
25 voltmeter 32. Optionally, a microprocessor 34 may be 

incorporated in the apparatus for summing the signals from the 
scaling amplifiers 24, 26 and/or further processing those 

30 

signals for diagnostic purposes, all in accordance with 
techniques known to those of ordinary skill in the art. If the 
35 microprocessor 34 is used to sum the signals from the scaling 
amplifiers, the summing amplifier 28 may be eliminated. 
The present invention is for an improved method for 

40 

calibrating the signals from the rib cage and abdomen such that 
their sum produces a signal indicative of tidal volume. As is 
45 noted above, while various calibration techniques are known to 
those of ordinary skill in the art, all the known techniques 
possess drawbacks. 
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5 Turning to the calibration technique of the invention, from 
Equation A it. is known that: 

V * (KxRC)+(LxAB) 

10 

where V is total respiration volume or tidal volume, RC is the 
rib cage contribution to respiration volume as measured by the 
15 rib cage conductor and AB is the abdominal contribution to 
respiration volume as measured by the abdominal conductor. K 
and L are calibration factors for the rib cage end "abdomen, 

20 

respectively. Another way of expressing this relationship is: 

V » Mxr (ZxRC)+AB3 [EQUATION D] 

25 where MxZ is equal to K and M is equal to L. It will be 

apparent that Equation D separates the calibration components 
into a preportionality factor Z and a scaling factor M . Using 

30 

this approach, calibration may be viewed as a two step process. 
The first step is the determination of the correct 
35 proportionality factor Z satisfying Equation D, such that 

V - (ZXRCJ+AB [EQUATION E] 

In other words, the proportionality factor Z defines the 

40 

correct relative contributions of the rib cage (RC) and abdomen 
(AB) to tidal volume (V) . Classically, determination of the 
45 proportionality factor Z is determined by an isovolume 

calibration technique in which the subject breathes against a 
closed airway i.e. with no volume movement at the mouth, whereby 

50 

V=»0. In other words, during an isovolume manuver, the only 

movement of volume is between the rib cage and abdomen 
55 compartments, i.e. pendelluft, since no air escapes through the 

mouth. Under these conditions, Equation D becomes: 
go O-(ZxRC) + AB [EQUATION F] 

or 



5 So, 'BJ^racording the readings from the rib cage and abdomen 
conductors during the isovolume manuver, the proportionality 
factor Z can be determined from Equation G. Once Z is 

10 

determined, the quantity (ZxRC)+AB can be calculated for any 
point in time from the recorded values of the rib cage and 
15 abdomen conductors. Since we know from Equation E that 
(ZxRC)+AB is always proportional to tidal volume V, a 
determination of that quantity provides a valuable diagnostic 

20 

tool. For example, as those of ordinary skill in the art will 
appreciate, from this quantity obstructive and central apneas 
25 can be diagnosed, RC as a percent of tidal volume V can be 

calculated, and increases and decreases in relative tidal volume 
V can be assessed* 

30 

The difficulty with this approach is that an isovolume 
manuver requires breathing against a closed airway, which is not 
35 always practical, as in neo-natal and critical care 
applications . 

40 
45 
50 
55 
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5 if, on the other hand, tha proportionality factor z and hence 
the quantity(ZxRC)+AB is to be determined without an isovolume 

iQ manuver, it would appear from Equation E that a measurement of 
tidal volume must be taken, otherwise there will be a single 
equation with two unknowns, namely, proportionality factor Z and 

15 tidal volume V. 

In accordance with the present invention, this problem is 
solved as follows. The rib cage (RC) and abdomen (AB) signals 
are recorded for a large number of breaths during an initial 
baseline period. For example, 250 breaths may be measured 
25- during quiet breathing over a 10 minute interval. Actually, 

breathing during the baseline period need not be quiet, as long 
3Q as it is steady state. For example, the baseline could be 
derived from breaths recorded during 10 minutes of exercise. 
During this baseline period, the uncalibrated signals from the 
35 rib cage (RC) and abdomen (AB) conductors are recorded. 

Referring to FIG. 3, for each of these signals, two values or 
4Q breath "deltas" are calculated for each breath, ( J being 
the difference between the signal at the beginning and end of 
inspiration, the other ( 2 ) being the difference between the 
45 beginning and end of expiration. These delta values are then 

totaled separately for each signal. Assuming 250 breaths during 
_ q baseline, and since there are two delta values for each breath, 
the total for each signal will be computed by adding 500 delta 
values. Where the delta values described above are preferred, 
55 it should be appreciated that the delta values are employed to 
provide a parameter indicative of the relative amplitude for 
each breath of the uncalibrated rib cage and abdomen signals 
taken during baseline. Accordingly, as used herein, the term 
delta values means any parameter of the uncalibrated rib cage 
65 and abdomen signals which provides this information. 
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5 if actual tidal volu*-* were also recorded during baseline, 

as by- spirometry, and the deltas for tidal volume totaled and 
the mean determined, the following relationship would apply: 

10 

Mean SP ■ « Mean RC + Mean AB [EQUATION H] 

Where SP is actual tidal volume as determined, e«g. by 
15 spirometry. Since the mean values in Equation H are derived 
from uncalibrated signals , calibration factors are required if 
the lefthand side in Equation H is to equal the righthand side. * 

20 

The standard deviations (SD) of the mean values for the 
25 tidal volume (SP) rib cage (RC) and abdomen (AB) signals can be 
calculated. Herein, these will be expressed, respectively, as 
SD (SP), SD (RC), and SD (AB) . If tidal volume (SP) is constant 

30 

for all breaths recorded during baseline, Equation H can still 
be computed and the standard deviation of tidal volume, SD (V) , 
35 is 0. 

This situation is analogous to Equation F which, as 
explained above, applies to the isovolume manuver where V=»0. In 

40 

particular, by considering the standard deviation of a constant 
tidal volume SD (V), which is also 0, the pendulluft occurring 
45 during normal breathing creates a situation analogous to the 
isovolume situation of Equation F, and taking the variance (Var) 
of both sides, yields: 

50 

Var (ZxRC) ■ Var (AB) [EQUATION J] 

Where variance (Var) is equal to (SD) 2 . Equation J may be 
55 expressed as 

Z 2 xVar (RC) = Var (AB) [EQUATION K] 

By taking the square root of both sides 

60 

ZXSD (RC) = SD (AB) [EQUATION L] 

Which yields 

65 Z = SD (AB) / SD (RC) [EQUATION M] 
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5 should be appreciated that the standard deviations of the 

means of the rib cage and abdomen signals is indicative of the 
average variability of those signals. Accordingly, 1 any 

10 

analysis that provides an indication of average variability of 
those signals may be used in lieu of computing standard 
15 deviations. 
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5 put into other words, it will be apparent from Equation M that 
the proportionality constant Z for solution of Equation E can be 
calculated, assuming constant tidal volume breathing during 

10 

baseline, from the ratio of the standard deviations of A6 and RC 
in a manner analogous to the isovolume manuver, but without the 
15 requirement that the airway be blocked. 

Of course, taking spirometry readings with a spirometer 
during baseline defeats an objective of the calibration 

20 

technique in accordance with the invention, namely, to calibrate 
non-invasively . 

25 However, baseline constant tidal volume assumption required to 
practice the calibration technique of the invention without 
invasively recording actual tidal volumes can be satisfied by 

30 

removing wild points from the 500 delta values of AB and RC 
computed during baseline as by eliminating values outside of 1.5 
35 standard deviations of the uncalibrated sum of the RC and AB 
components. With those delta, values excluded, the remaining 
delta values for RC and AB are separately totaled, the means 

40 

determined, and the standard deviations for the means 
calculated. The proportionality factor Z can then be calculated 
45 from Equation M. 

Once the proportionality factor Z is known, the quantity 
(ZxRC)+ AB will always be proportional to actual tidal volume. 

50 

See Equation E. This quantity can be continuously monitored on 
a real time basis from the real time rib cage and abdomen 
55 signals generated by the apparatus of FIG. 2. Preferably, this 
quantity is expressed as a percent of (ZxMean RC)+Mean AB where 
Mean RC and Mean AB are, respectively, the means of the totals 

60 

of the uncalibrated RC and AB delta values generated at 
baseline, but exluding the wild points. This is sufficient for 
65 the bulk of diagnostic work, such as detecting obstructive and 
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3 central apneas, hypoapeneas, and variations in tidal volume, the 
latter being important as a diagnostic tool in a wide variety of 
disorders . 
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Referring to FIG. 2, assuming the scaling amplifiers have 
initially been set to unity gain, once the proportionality 
factor z is determined, the scaling amplifier 24 for the rib 
cage is adjusted to Z. From Equation E, we then know that the 
sum of the signals from the scaling amplifiers at the output of 
15 the summing amplifier 28 will be preportional to tidal volume. 

This completes the calibration procedure. 
2Q While for most purposes a guantatitive determination of 
tidal volume is not necessary, guantitive calibration can be 
achieved once the proportionality factor Z is known, in 
25 particular, referring to Equation D, by taking a single 

measurement of actual tidal volume V, the scaling factor M can 
be calculated, since ZxRC and AB are available, respectively, at 
the output of the rib cage arid abdomen scaling amplifiers, in 
other words, there is then only a single unknown in Equation D, 
35 the scaling factor M, which is calculated as: 

M = V/[(ZxRC)+AB]. [EQUATION N] 

One simple way to take a measurement of actual tidal volume is 

40 

to simply have the subject inhale a known quantity of air, as 
from a syringe. Before the subject inhales, this quantity is 

45 input to the apparatus of FIG. 2 at the "spirometer", input. 
The microprocessor 34 can then perform the calculation of 

^ Equation M from the values of Z, RC and AB at the end of 
inspiration from the syringe. 

The scaling factor M may be set in the FIG. 2 apparatus by 

55 multiplying the gain of the summing amplifier 28 by scaling 
factor M, whereupon the output of the summing amplifier will be 
a semi-quantatitive indication of tidal volume. The term 

60 

semi-quantatitive is used because it has been determined that 
tidal volume computed in this fashion is ± 10% of actual tidal 
65 volume as 
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3 determined by spirometry. 

Desirably, the accuracy of the proportionality factor Z is 
iq continuously monitored by repeatedly recalculating the 

proportionality factor at five minute intervals and printing out 

the resulting value. After the initial calibration procedure, 
15 the recalculated value for the proportionality factor Z during 

each subsequent five minute interval should be 1.0. In the 
2Q event of a change in the position of the patient or other 

condition that varies the proportionality factor Z, that will 

exhibit itself as a reading of z above or below 1.0. If the- 
25 variations are too large, i.e. if z is less than about .7 or 

greater than about 1.3, the calibration routine can be rerun. 
3q Suitable visual or audio alarms can be incoporated in the 

apparatus of FIG. 2 to indicate such excessive variations. 



40 



45 



50 



55 



60 



WO 87/04332 



PCT/US87/00217 



-18- 

5 j^noted above, the calibration technique in accordance with 
the present invention is currently utilized in connection with 
an apparatus for monitoring respiration sold by Nims, Inc., 
Miami Beach, Florida, under their model designation Respigraph 
\ The Respigraph incorporates a microprocessor. 

15 Preferably, the microprocessor is programmed to carry out the 
calibration technique in accordance with the present invention. 

20 A program listin 9 tor a suitable program for carrying out the 
calibration function appears below in Table A. The program is 
in assembly language for a Z80 microprocessor as manufactured by 

25 Zylog, Inc. 
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. I DENT GDC 
. INSERT SP80.ASM 
. INSERT FPMAC . SRC 
. INSERT CALCOM.ASM 
.INSERT CKCOM.ASM 
. INSERT UTC'COM 
. INSERT 3ELC0M .ASH 
. INSERT VARCOM.CCC 
. INSERT GL3C0M 
.INSERT 3YC0M 
. INSERT PKCOM 

. INSERT EFFCOM 
. INSERT SIOCCM 
. INSERT HRDCOM 

. EXTERN 'v'ALCOL ,','ALFND , TMCN , SAVECL , LPON , LPOFF .^UOFF 

.EXTERN FIND 

. EXTERN DERIV5 .PRTIME 

.EXTERN SAUECL 

.EXTERN UAL 



■ENTRY GDC1 ,QDC2.L0DC1 
. ENTRY BASE 

. EXTERN RTQN.RTQF.RTICK 

.EXTERN QOCAL,SYST,ANSOUT 



. B<TERN SRTINIT.SRTPLOT.EXEC 
•EXTERN LABEL 
. EXTERN TDHE 

.EXTERN PTRUPO , GETP2 , ROUND 



FINIT 



QDC: 
i 

;THESE JUMPS MUST BE THE 1ST THINGS: 

JMP QDC1 

JMP QDC2 

JMP UAL 

JMP LQDC1 

JMP ISGV1 

JMP BASE 
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ISQUli 



BASE: 



XRA A 

STA BSEFLG 

MVI A.QFFH 
STA . ISOFLG 

MUI A,0 

STA WTFLG 

JMP QOCCOL 



mi A.OFFH 

STA BSEFLG 

MVI A ,OFFH 

STA WTFLG 

XRA A 

STA ISOFLG 

JMP Q0C3C0L 



QDC l : 



LQDC1 



XRA 
STA 
MUI 
STA 
STA 
JMP 



XRA 
STA 
MVI 
STA 
XRA 
STA 
JMP 



A 

8SEFLG 

A,0 

WTFLG 

ISOFLG 

QOCCOL 



A 

BSEFLG 
A,OFFH 
WTFLG 
A 

ISOFLG 
QOCCOL 



QOCCOL: 



XRA 
CMA 
STA 
XRA 
STA 



QUALFLG i QUALITATIVE 
A 

BSEFLG 



WO 87/04332 
FILE! QDCASM 



21 



PCT/US87/00217 



! 

QDC3C0L : 



START THE DATA COLLECTION 



XRA A 
STA MALFLG 

STA P0S1 ;NOT IN SINGLE POSITION 
LXI H,0 

SHLO NBERS ;t» VALUES IN -VALUES 4 YUALUES 

i 

LXI H,0 
SHLO TIME 
LKI H,0 
SHLD MAXCNT 

• 

LDA BSEFLG 

IF <.A,IS,ZEROi ;NGT COLLECTING 8ASELINE . 
jINIT THE TRANSDUCESERS 
ml A, 10 
CALL SYST 

•THIS RESETS CALS TO 1 AND 1 30 AFTER THIS IS EXECUTED, ALL OTHER 
; CALIBRATIONS ARE NULLIFIED 



LOOP 



XRA 


A 


STA 


CALTYPE 


STA 


LASTCAL ; RIGHT NOW UNCAL8RATED 


STA 


GOODBASE ;N0 LONGER HAS 


CALL 


LPON 


FPRN 


RES ,RESPT 


CALL 


LPOFF 


ENDIF 




CALL 


SETUP 


CALL 

i 


SRTINIT ; IN IT RT. PLOT 


CALL 


SRTPLOT 


LDA 


BFLG 


IF <.A, 


IS.NZERO) 


: # DELTAS AND ELAPSED TIME 


MVI 


A, 3 


STA 


OPCODE 


CALL 


ANSOUT 


ELSE 




LDA 


SECFLG 


IF <.A, 


IS.NZERO) 


mi 


A,3 


STA 


OPCODE 
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CALL ANSOUT 
LDA FLSHFL6 
OCR A 
STA FLSHFLG 
IF <.A,IS,ZERO> 

LXI .H.TOPLFT+53 
LDA BSEFLG 
IF <. A, IS, ZERO) 
LDA I50FLG 
IF (.A. IS, ZERO) 
LXI O.CLBMESS 
ELSE 
LXI D.ISOtlESS 

END IF 
ELSE 

LXI D.BSEMESS 
ENDIF 

CALL LABEL 
MVI A, 2 
STA FLSHFLG 



ELSE 



LXI H,T0PLFTt53 
LXI D, BLANKS 

CALL LABEL 



ENDIF 

Er* - 
ENDIF 

LDA GFLG ;= OFFH THEN CLOCK 

jSTILL ON 
EXITIF< .A,IS,ZERO) 
ENDLOOP 



NOW DATA COLLECTION PHASE HAS BEEN COMPLETED 
MAXCNT HOLDS THE NUMBER OF MAX POINTS COLLECTED 

IRSRV USES CHANNEL SPBUFF FOR DETECTIONS 
SPMAX HOLDS THE DELTAS FROM SPBUFF 
RCMAX HOLDS THE DELTAS FROM RCBUFF 
ABMAX HOLDS THE DELTAS FROM ABBUFF 

CLEAR THE SCREEN 

CALL TVON 
CALL TVOFF 



LDA BSEFLG 
IF ( .A, IS, ZERO) 
LDA ISOFLG 

IF <.A,IS,NZERO) ;THIS WAS AN ISGVOLUME MANUEVER 

JMP ISOCALC 
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ENDIF 
ENDIF 

LHLD MAXCNT 

LXI 0,3 

ORA A 

DS8C 0 

IF <PSU, IS, CARRY) 

LDA WTFLG 
IF <. A, IS, ZERO) 
ffJI A,l 

STA OPCODE jNOT ENOUGH DELTAS 
CALL ANSOUT 

mi a,o 

STA OPCODE i PR INT MENU ON SCREEN 
CALL ANSOUT 
ENDIF 
STC 
RET 
ENDIF 



jTHERE ARE 2*<MAXCNT -1) FLOATIND PT VALUES IN X^LUES 
i AND YVALUES FOR THIS RUN 
LHLD MAXCNT 

SHLD N8ERS jTOTAL » STORED 
CALL QDCAL 
i RETURN WITH CARRY SET IF BAO FOR ANY REASON 



ORA A 
MVI A,0 

STA OPCODE {PRINT MENU ON SCREEN 
CALL ANSOUT 
RET 



ISOCALC: 

LDA URAPF 

IF ( .A,IS,N2ER0) jTHEN 400 PTS AND RIPT PTS TO 1ST 

j DATA PT 

LXI H.400 
SHLD CNT 
LHLD RIPT 
SHLD ROPT 
ELSE 

iTIME HOLDS » OF DATA POINTS AND 1ST PT AT OFFSET = 0 
LHLD TIME 
SHLD CNT 
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L00P1 : 



IX I 


H,0 


SHLO 


ROPT 


ENOIF 




FFOA ' 


0 


FSTA 


XSUM 


FSTA 


YSUM 


FSTA 


XXSUM 


FSTA 


XYSUM 


FSTA 


XCNT 


CALL 


GETX 


CALL 


GETY 


CALL 


XYPAIR 


LHLD 


ROPT 


LXI 


B.RSIZE 


CALL 


PTRUPD 


SHLD 


ROPT 


DSKZ 


CNT.L00P1 


CALL 


LSTFIT 


FLDA 


SLOPE 


FABS 




FSTA 


ABCAL 


FFOA 


1 


FSTA 


RCCAL 


IT ON SCREEN 


MU1 


A 1 4 


STA 


OPCODE 


CALL 


ANSOUT 


CALL 


SAVECL 


IF <PSW 


,IS,NCARRY) 


MVI 


A, 3 ;I 


STA 


CALTYPE 


mi 


A.OFFH 


STA 


LASTCAL ;G 


XRA 


A 


cm 




STA 


QUALFL6 


CALL 


LP ON 


CALL 


PRTIME 


FPRN 


1SV.ISVPT 


CALL 


LPOFF 


ORA 


A 


END IF 





;THEN GOOD VALUES 
;ISOV QUAL 



{GOOD CALIBRATION 



RET 

5 

t 
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ISV: .ASCIS •' 

QUALITATIVE ISOVOLUME CALIBRATION 
RCCAL = 303 
ABCAL = 3D3 

ISVPT: .WORD RCCAL, ABCAL 



GETX: 



LXI 

LDEO 

DAD 

MOV 

I NX 

MOV 

SDED 

FILA 

FLOT 

FSTA 

RET 



H.ABBUFF 

ROPT 

D 

E,M 

H 

D,M 

XVAL 

XVAL 



GETY ! 



LXI 

LDED 

DAD 

MOV 

I NX 

MOV 

SDED 

FILA 

FLOT 

FSTA 

RET 



H.RCBUFF 
ROPT 

D 

E,M 

H- 

D,M 

YVAL 

YVAL 

YVAL 



XYPAIR: 



FLDA 


XVAL 


FLCfB 


XSUM 


FAOO 




FSTA 


XSUM 


FLDB 


YVAL 


FLOA 


YSUM 


FAOO 




FSTA 


YSUM 


FLOB 


XVAL 


FLDA 


XVAL 


FMUL 




FLDB 


XXSUM 


FAOD 




FSTA 


XXSUM 
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FLD8 XVAL 

FLDA YVAL 
FMUL 

FLDB XYSUM 
FA DO 

FSTA XYSUM 

FLDB XCNT 
FFOA 1 
FAOD 

FSTA XCNT 
RET 

LSTFIT: 

FLDA XSUM 
FLDB XSUM 
FMUL 
FPSH 



FLDA XCNT 

FLDB XXSUM 

FMUL 

FSWP 

FPOP 

FSUB 
FPSH 

FLDA XSUM 
FLDB YSUM 
FMUL 
FPSH 

* 

FLDA XCNT 

FLDB XYSUM 

FMUL 

FSWP 

FPOP 

FSUB 
FSWP 
FPOP 
FDIV 

FSTA SLOPE 
RET 



QDC2: 

;THIS CAN BE DONE AFTER EITHER A Q0C1 (CALTYPE = 1/2) 
; OR AFTER AN ISOVOLUME < CALTYPE = 3/4) 

• 

•IF CALTYPE = 0, THEN HAVEN'T CALIBRATED THE RC TO AB 
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LDA CALTYPE 
IF (.A, IS, ZERO) 

MVI A, 2 

STA OPCODE 

CALL ANSOUT 
RET 
END IF 



iNOU COLLECT DATA AS IF VALIDATING 
CALL VALCOL 



LHLD 
LXI 
ORA 
DSBC 



MAXCNT 
0,2 
A 
D 



IF (PSW, IS, CARRY) 



MVI 
STA 
CALL 
WI 
STA 
CALL 
RET 
ENDIF 



A,i 

OPCODE 
ANSOUT 
A,0 

OPCODE 
ANSOUT 



i'THEN LESS THAN 2 BREATHS 



i PRINT MENU ON SCREEN 



NOW FIND THE RATIOS OF SUM/KNOWN VOLUME 

CALL VALFND 
RETURNS WITH AR£G=* OF VALUES USED IN SUMMATION 

FSTA FCNT 



NOW XVALUES HOLDS ALL OF THE RATIOS, XSUM 'HOLDS SUM OF ALL XVALUES 
XXSUM HOLS SUM OF SQUARES OF ALL XVAUES, AND FCNT HOLS THE NUM8ER 
OF VALUES USED IN THE SUMPTION 

COMPUTE THE MEAN RATIO 
FLDB XSUM 
FLDA FCNT 
FDIV 

FSTA MEAN 

i'.NOW KNOW THAT SUM NEEDS TO 8E MULTIPLIED 8Y 1/MEAM, SO 
! RCCAL - CURRENT RC GAIN (RCGNF) # 1/MEAN ' 
! ABCAL » CURRENT AB GAIN (ABGNF) » 1/MEAN 



FFOA 


1 


FATB 




FLDA 


MEAN 


FDIV 




FLDB 


RCGNF 


FMUL 




FSTA 


RCCAL 


FFOA 


1 
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FAT8 
FLDA 
FDIV 
FLDB 
FMUL 
FSTA 



MEAN 

ABGNF 

ABCAL 



MVI A, 4 
STA OPCODE 
CALL ANSOUT 

CALL SAVECL 
IF (PSW.IS.NCARRY) 
CALL LPON 
CALL PRTIME 
FPRN VL2.VL2PT 
LPOFF 
H,0 



iTHEN GOOD 



CALL 
LXI 

SHLO NBERS 
XRA A 

STA QUALFLG ;NOT QUALITATIVE, IT IS QUANTITATIVE CAL 



LDA CALTYPE 
IF <.A,LT,3) 
MVI A, 2 



STA 
MVI 
STA 
ELSE 

mvi 

STA 

mi 

STA 
ENDIF 



CALTYPE 
A,OFFH 
LAST CAL 



;THEN QOC METHOD 
;NOW ITS QDC QUANT 



:QOOD CAL 
(ISOVOLUME METHOD 

A 4 

CALTYPE jNOT ITS ISGV QUANT 
A.OFFH 

LASTCAL ;GOOD CAL 



5MULT THE CURRENT VALUE OF VTLIMIT BY 1/MEAN SO WILL 8E SCALED OK 
; NOTE: VTLIMIT IS INITIALLY SET UHEN DOING QDC1 

{THIS IS FOR INITIALLI2ING THE HYPOPNEA VOLUME FOR THE 1ST 5 MINUTES! 
FLDB VTLIMIT 
FLDA MEAN 
FDIV 

FSTA VTLIMIT 

)*»»**»*»*****»«♦*##*»»«»#*»*»»»*»««***#*#*»»»♦»**»**♦»*******»* 
5 PROBLEM! UHAT IF DID ISOVOLUME 8EF0RE HERE, THEN VTLIMIT WAS NOT 
; EVER INITIALIZED ???? 

;»*♦»**»♦***»***»♦«**♦************»***♦»»«*»»****♦*****»»***»*** 
ENDIF 
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FILE: QDC.ASW 



>mvi A,0 

jSTA OPCODE {PRINT MENU ON SCREEN 
(CALL ANSOUT 
RET 



VL2: .ASCIS ' 

QUANTITATIVE CAL FACTORS 
RC CAL = 3D3 
AB CAL = 3D3 

YL2PT I 

.WORD RCCAL.ABCAL 



;NOTE: MCA VALUES OF RCINVOL, ABINVOL, AND SMVOL UILL 

! BE STORED BEGINNING AT ADDR = XVALUES 

! XVALUES + YVALUES = 1400 BYTES 

I I UILL STORE UP TO 1S00 8YTES 

! = > 500 DELTAS = 250 INSP DELTAS * 6 BYT-5 

! » 1500 BYTES 

! MCAADDR ■ ADDRESS OF NEXT TO STORE 



SETUP: 

jTHIS IS ONLY USED WHEN CALIBRATING QOC METHOD 
j SETUP THE CLOCK 

} SAMPLE AT 20 PTS PER SEC 

PUSH H 
i 1NIT REAL-TIME BUFFER POINTERS 
; MODIFIED 09-24-81 
iFOR PEAK FLOWS : 

LXI H.DLTIMES 

SHLD TMPTR 

LXI H,0 

SHLD TTPOS 

SHLD TTNEG 

SHLD TEMPK 
• SHLD TTPTM 

SHLD TTNTM 

SHLD PKTIME 

LXI H,0 

SHLD DLTIPTR j OFFSET TO STORE NEXT DELTA 

XRA A 

STA DLTFULL j DELTA RING BUFFER NOT FULL YET 
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LTTR 






STA 


FSTMN 

• 




LXI 


H.XVALUES ;WILL USE THIS BUFFER TO 


<» 


SHLD 


MCAADDR ;MCA VALUES 




LXI 


H,0 




SHLD 


SPTEMP 




XRA 


A 




STA 


FRSTMX ;1ST MAX NOT FOUND 




LXI 


H.TMBUFF 




mi 


M,0 


* 


I NX 


H 




MVI 


M , 0 




MVI 


A, 2 




STA 


FLSHFLG 




XRA 


A 




STA 


WRAPF 




STA 


SECFLG 




CMA 






STA 


BFLG 




LXI 






SHLD 


ELPTM 




SHLD 


ELPMIN 




SHLD 


SECONDS 




LXI 


H,0 




SHLD 


RIPT 




SHLD 


ROPT 




SHLD 


ODRIPT 




LXI 


H,0 




9HLD 


TIME 





; 1NIT VARIABLES FOR PICKING MAX AND MINS 



XRA 


A 




STA 


MAXB 




STA 


PTF 


; START BY LOOKING FOR 


STA 


SLOPEF 


jINIT MAX NOT FOUND 


XRA 


A 




STA 


INTFLG 




CMA 






STA 


GFLG 


{INDICATES CLOCK ON 


LXI 


H,IRSRV 




CALL 


RTON 


;TURN ON CLOCK 


POP 


H 
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RET 



CLKSTP s 

iSTOP DATA COLLECTION 
{TURN CLOCK OFF 

MODIFIED 09-24-81 

CALL RTOF {TURN OFF CLOCK 

XRA A 

STA GFLG 

POP D 

POP H 

POP B 

POP PSU 
RET 



IRSRV: 
I 

; INTERRUPT SERVICE ROUTINE 



PUSH 


PSW 


PUSH 


B 


PUSH 


H 


PUSH 


D 


CALL 


RTICK 


LDA 


INTFLG 


IF <.A 


IS.NZERO) 


NOP 




ELSE 




CMA 




STA 


INTFLG 


ENDIF 




LHLD 


ELPTM 


I NX 


H 


SHLD 


ELPTM 


LDED 


SFREQ 


CALL 


TDHE 


IF <PSU 


, IS, ZERO) 


LXI 


H,0 


SHLD 


ELPTM 


LHLD 


SECONDS 


INX 


H 


SHLD 


SECONDS 


XRA 


A 



{SERVICE INTERRUPT 

{THEN INTERRUPT INTERRUPT 

{SET FLAG 



{THEN 1 SECON 
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CMA 

STA SECFLG 
LXI D,60 
CALL TDHE 

IF <PSW,IS,ZERO> ;THEN 1 MINUTE 

LXI H,0 

SHLD SECONDS 
LHLD ELPMIN 
INX H 

SHLD ELPMIN 
;IF COLLECING BASELINE: THEN STOP AT 5 MINUTES 
LDA BSEFLG 
IF < .A.IS.NZERO) 
LXI D,5 
LHLD ELPMIN 
CALL TDHE 
IF <PSU,IS,2ERO) 

JMP CLKSTP 

END IF 
ELSE 



WHILE CALIBRATING: IF REACH 10 MINUTES BEFORE 500 DELTAS, 
THEN AUTOMATICALLY STOP 

LXI D , 1 0 

LHLD ELPMIN 
CALL TDHE 

IF (PSW.IS.ZERO) ; HAVE REACHED 10 MINUTES 



JMP CLKSTP 
END IF 
ENDIF 



ENDIF 
ENDIF 



ALWAYS COLLECT RC INTO RCBUFF, 
; AB INTO ABBUFF, 

NOTE: ALWAYS USES DATA IN SPBUFF FOR DETECTION 

j RIB CAGE VALUE 



LHLD 


RCAD 


MOV 


A , H 


cm 




AN I 


OFH 


MOV 


H , A 


MOW 


A,L 


CMA 




MOV 


L,A 


SHLD 


DATA 


SHLD 


RCV 
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LHLO RIPT 
LXI D.RCBUFF 
CALL PRNG 



LHLD A BAD {ABDOMEN UALUE 

MOV A,H 

CMA 

AN I OFH 

MOV H,A 

MOV A,L 
CMA 

MOV L,A 

SHLD DATA 

SHLD ABV 

LHLD RIPT 

LXI D.ABBUFF 

CALL PRNG 

; 

; SUM INTO SPBUFF TO PICK OFF BREATHS 

LHLD RCV 

LDED ABV 

DAD D 

SHLD DATA 

! 

LHLD RIPT 

LXI D, SPBUFF 

CALL PRNG 



i INCREMENT RIPT FOR NEXT TIME 



LHLD 


TIME 


I NX 


H 


SHLD 


TIME 


SHLD 


DATA 


LHLD 


RIPT 


LXI 


D.TMBUFF 


CALL 


PRNG 


LHLD 


RIPT 


SHLD 


ODRIPT 


LXI 


B.RSIZE 


LHLD 


RIPT 


CALL 


PTRUPD 


SHLD 


RIPT 
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LDA URAPF 

IF <.A,IS,ZERO> | SEE IF RIPT RESET TO 0 
LXI D,0 
LHLD RIPT 
CALL TDHE 
IF <PSU,IS,ZERO> 
MVI A.OFFH 
STA URAPF 
ENDIF 
END IF 



LDA BSEFLG 
IF <.A,IS,2ER0) 
LDA ISOFLG 
IF CA.IS.NZERO) 

JMP EXIT 
ENDIF 
ENDIF 



iTHEN NO PEAK DETECTION 



LDA SLOPEF 

IF C.A,IS,2ER0) j THEN INIT POINTS NOT COLLE ITED YET 
CALL GETP2 ; RETURNS C IF < 5 PTS IN B'.-FER 

{THIS JUST MAKES S.RE EXTRA 
; POINTS IN BUFFER 
JC EXIT fNOT ENOUGH POINTS IN BUFFER 

NOTE: GETP2 INITIALLIZES ROPT TO LINE UP WITH DERIY STUFF 
SKIP 1ST 3 POINTS IN BUFFER 

LXI H,4 ; OFFSET POINTING TO POINT 4 

SHLD ROPT 



MVI A.OFFH 
STA SLOPEF 
ENDIF 



CALL DERIV5 



MUST ADD THE COOE TO DO PEAK DETECTION 
HERE TOO ! ! 

LDA FRSTMX 

IF < .A, IS, ZERO) ;HAVE NOT FOUND 1ST MAX POINT 
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LHLD ROPT 
CALL FIND 

IF <PSW,IS,NCARRY) ;N0T A MAX 
JMP OUTER 

ENDIF 
;HAVE 1ST MAX 

LDA FRSTMX 
CMA 

STA FRSTMX 

SHLD SPTEMP ;TEMP CRITICAL POINT 
SDED RCTEMP 
SBCD ABTEMP 
LDA RTF 

cma 

STA PTF ;NOU LOOK FOR MIN 
JMP OUTER 
ENDIF 



LHLD ROPT 

CALL FIND 

IF (PSW, IS, CARRY) 
SHLD SPVALUE 
SDED RCVALUE 
SBCD ABVALUE 



LDA 


i RTF 




IF 


< .A, IS, ZERO) ;THIS IS A MAX POINT 




LHLD 


RCVALUE 




SHLD 


RCMXVAL 




LHLD 


ABVALUE 


i 


SHLD 


ABMXVAL 


1 


CALL 


OLTSTORE 


;ONLY STORE 


THE MCA VARIABLES AT A MAX < ONLY WANT FOR 


• 

• 


CALL 


MCASTORE 


• 

llF DLTIPTR 


= 1000, THEN EXIT 




LHLD 


DLTIPTR 




LXI 


0,1000 




CALL 


TDHE 




IF <PSU 


, IS, ZERO) 




mvi 


A , OFFH 




STA 


DLTFULL 




LXI 


H,0 




SHLD 


DLTIPTR 




JMP 


CLKSTP 




ENDIF 




ELSE 


iTHIS IS MIN POINT 




LHLD 


RCVALUE 




SHLD 


RCMNVAL 
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LHLD ABVALUE 

SHLD ABMNVAL 

LDA FSTMN 

IF < .A,IS,NZERO> ,'THIS IS THE. 1ST MIN, DO NOT STORE 

LDA FSTMN 
CMA 

STA FSTMN 

ELSE i STORE THE DELTA 

CALL OLTSTORE 



;IF DLTIPTR = 1000, THEN EXIT 
LHLD DLTIPTR 
LXI 0 ,1000 
CALL TDHE 
IF (PSW,IS,ZERO) 
MY I A.OFFH 
STA DLTFULL 
LXI H,0 

SHLD DLTIPTR 
JMP CLKSTP 
EN01F 
END IF 

END IF 
LDA PTF 

CMA 

STA PTF 
END! F 

t 
I 

OUTER; 

j UPDATE ROPT FOR NEXT TIME 
• 

LHLD ROPT 
LXI B.RSI2E 
CALL PTRUPD 
SHLD ROPT 



EXIT: 

; LOOK FOR ENTER 

j 

call csts 
if •:psw,is,nzero) 
call cinp 

if <.a,eq,ooh) ;if enter 

jmp clkstp 

END IF 
END I F 

• 

XRA A 



i 
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STA INTFLG 

POP 0 

POP H 

POP 8 

POP PSU 
RET 



MCASTORE: 

LHLD MCAAOOR 

LOED RCINVOL 

MOU M.E 

I NX H 

MOV M,D 

INX H 

LDEO ABINVOL 

MOV M,E 

INX H 

MOV M,0 

INX H 

LDEO SMVOl 

MOV M , E 

I NX H 

MOV M,0 

INX H 

SHLD MCAAOOR 
RET 



OLTSTORE : 

LHLD 
LDED 
ORA 
DSBC 

;HL = 1NSP VOLUME 

SHLD 

LHLD 

LXI 

DAD 

LDED 

MOV 

INX 

MOV 



LHLD 
LDED 
ORA 
DSBC 
!NSP VOLUME 

SHLD 
LHLD 
LXI 



;HL 



RCMXVAL 

RCMNVAL 

A 

D 

jRC 

RCVALUE ; DELTA RC 
DLTIPTR j OFFSET TO STORE 
D , RCDATA 
D 

RCVALUE 

M,E 

H 

M,D 

A6MXVAL 

ABMNVAL 

A 

0 

;A6 

A8VALUE jDETLA AB 
DLTIPTR (OFFSET TO STORE 
D .ABDATA 
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DAD 0 

LDED ABVALUE 

MOV M . E 

I NX H' 

MOV M,0 

LHLD DLTIPTR 

I NX H 

INX H 

SHLD DLTIPTR 

LHLD MAXCNT 

INX H 

SHLD MAXCNT 

XRA A 
CMA 

STA BFLG 
i STORE THE TIMES 



LHLD 


TMPTR 


LXI 


D , DLTEND 


CALL 


TDHE 


IF <PSW 


,IS,2ERO) 


RET 




END IF 




LHLD 


TMPTR 


LDED 


CPTIME 


MOV 


M,E 


INX 


H 


MOV 


M , D 


INX 


H 


LDED 


PKVALUE 


MOV 


M,E 


INX 


H 


MOV 


M,D 


INX 


H 


SHLD 


TMPTR 




RET 



;NO ROOM 



;TIME OF THE CP 



;TIME OF THE PEAK FLOW 



;FOR NEXT TIME 



PRNGs 



HL HOLDS THE RING INPUT PTR 
DE HOLDS ADDR OF RING 



DAD 

LDED 

MOV 

INX 

MOV 

RET 



! 

ISOMESS: 
i 

BSEMESS: 



D 

DATA 

M,E 

H 

M,D 



.ASCIS 'Isovolume' 
.ASCIS 'Baseline' 
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CL8MESS: 


.ASCIS 


'Cal ibrating' 


8LANKS: 
* 


.ASCIS 




RES: .ASCIS 






AT 31:31 i3I 


Cal -factors reset to 1 


RESPT: .WORD 


HOUR MIN 


1 . SEC 


i 
1 

.EXTERN 


QDCQM 

LI 1/ w Ul ■ 




.LOC 


QDCQM 




i 

MAXB: .BLKB 


1 




RCCU: .BLKB 


2 




ABCU: .BLKB 


2 




SPCU: .BLKB 


2 




RCV : . BLKB 


2 




ABV: .BLKB 


2 




MEAN : .BLKB 


4 




FCNT: .BLK8 
■ 


4 




FRSTMX : 


.BLKB 


1 


RCMXVAL : 


.BLKB 


2 


RCMNYAL : 


.BLKB 


2 


A3MXUAL: 


.BLKB 


2 


ABMNVAL: 


.BLKB 


2 


FSTMN: .BLKB 


1 




FLSHFLG: 


.BLKB 


1 


ISOFLG: 


.BLKB 


1 


XCNT: .BLKB 


4 




QOCEXTRA: 


.BLKB 


9 j EXTRA 


.RELOC 






.END 


QDC 
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. I DENT QDCAL 
. INSERT FPMAC.SRC 
. INSERT SPBO 
FINIT 

.INSERT VTCCOM 
. INSERT CALCOM 
.INSERT OLQCOM 
. INSERT CKCOM 
.INSERT 'JARCOM.CCC 
.INSERT SYCOM 
. INSERT EFFCOM 
. INSERT SIOCOM 
. INSERT QD1 .ASM 
. INSERT PKCC'M 

• 

. EXTERN CMPSTATS 

* 

.ENTRY QDCAL 

. EXTERN TDHE , ROUND 

. EXTERN SAUECL,ANSOUT,LP0N,LP0FF 

.ENTRY V2 iSQ CAN LOOK AT MAP! 

.EXTERN PRTIME 

.EXTERN WTSAVECL 

l 

.EXTERN TVQN.TVQFF 

. 

! 

; 
i 

QDCAL: 

LDA BSEFLG 

IF < .A.IS.N2ER0) (BASELINE 

JMP BCOMP 
END IF 

i ENTERED HERE. THEN IN CALIBRATE MODE (NOT LONG TERM' 1 
CALL CMPSTATS 

RC ; CARRY SET IF BAD CAL 

: BECAUSE NOT ENOUGH DATA COLLECTED 



;ABCAL = RCSTD 7 ABSTD 

FLDB RCSTD 

FLDA ABSTD 
FDIV 

FSTA ABCAL 
jRCCAL = i 

FFOA 1 

FSTA RCCAL 



MUST STORE THE GAINS, BECAUSE RESCALE ASSUMES SO !! 
LDA UTFLG 



( 
I 
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IF (.A, IS, ZERO) :WANT TO PRINT ERROR CODE IF BAD 

CALL SAUECL 
ELSE 

CALL WTSAUECL 
END IF 

IF '.PSW, IS, CARRY) iTHEN BAD CALS 

mi A.O jNO GOOD CALIBRATION (WHEN BEGAN THIS RUN. RESET 
STA CALTYPE |CALS TO 1 AND 1 
MWI A.O 

STA LASTCAL (NO GOOD CAL 
STA GOODBASE :N0 GOOD BASELINE NOW 

i PRINT OUT WHAT THE BAD CALS WERE AND THE STDS' 
CALL BDCALS 
STC iBAD CALS 

RET ^ 



;NOU MUST USE THE MEAN VT MEASURED TO SCALE THE CAL FACTORS 
; MEAN VT = VTMX/10 

CALL RESCALE 

LDA WTFLG 
IF .A, IS ,ZERO) 

MWI A ,4 

STA OPCODE 

CALL ANSOUT j DISPLAY THE CALS ON THE MvM T OR 
END IF 

LDA WTFLG 

IF (.A, IS, ZERO) iWANT TO PRINT ERROR CODE IF SAO 

CALL SAVECL 
ELSE 

CALL WTSAUECL 

IF <PSU,IS,NCARRY) (THEN GOOD CALS 

CALL LPON 

CALL PRTIME 

FPRN VI ,VtPT 

CALL LPOFF 

CALL CMPMCA ; COMPUTE MCA/VT . 

CALL LPON 

FPRN V2,V2PT 

CALL PRTHYP 

CALL LPOFF 

MVI A,l ;QDC DUAL 

STA CALTYPE 

MVI A , OFFH 

STA LASTCAL 

ORA A iGOOO CALS 

ELSE 

MVI A,0 ;N0 GOOD CALIBRATION (WHEN BEGAN THIS RUN, RESET 

STA CALTYPE ;CALS TO 1 AND 1 
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FILE i QDCAL.ASM 



MVI A,0 

STA LA?TCAL ;N0 GOOD CAL 

STA GOODBASE ;N0 GOOD BASELINE NOW 

; PRINT Q'JT WHAT THE 3A0 CALS MERE AND THE STD5 



CALL 
FFOA 
F5TA 
FSTA 

call 

CALL 

FPRN 

CALL 

LXI 

CALL 

CALL 



80CALS 
1 

RCCAL 

A8CAL 

CMPMCA 

LPON 

V2 .V2PT 

PRThYP 

H.STmP.S 

TXTYP 

LPOFF . 



i COMPUTE MCAA'T 



STC 
ENDIF 
RET 



i BAD CALS 



SDCALSi 



CALL LPON 

LXI H. STARS 

CALL TXTYP 

CALL PRTIHE 

FPRN BADC.BAOPT 

CALL LPOFF 



RET 



PRTHYP : 



FPRN HYP , HYPPT 
LDA QUALFLG 
IF < .A,IS,NZERO> 

LXI H.RHYP 
ELSE ; QUANT 

LXI H.AHYP 
ET-JDI F 
CALL TXTYP 
RET 



jQUAL 



BCQMP: 

;THIS COMPUTE MEAN VT AND IF I SOU WAS QUALITATIVE. CALT"?E = 3, 
; THEN IT SCALES IT TO YTMAX'IO 
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CALL CMP STAT 5 

IF -IPSW.IS.ZERO) . ;N0T ENOUGH OtfTA 

XRA A 

STA G0C03A3E [NOT A GOOD EASEL INE 

STC 
RET 
END IF 

t 

PESCALz i THIS WILL C:->-<NGE vrtLvE OF RCCAL AND A8CAL 

'-DA calt-fe 

IF •• .A.EQ.3) i ISC' J QUAL 

Isave cals 

snow must use the mean m? measured to scale the ^al fa~"rs so 

! MEAN VT = UTMX/iO ' 
CALL WT3A^=CL> 

i V. ... 

i 
! 

IF (PSU, IS.NChRR v > {THEN GOOD CALS 

, «V: A, 5 ! ISC -QUAL WITH BASELINE' " ' 
V^JTA CALTVP5 
MVI A.OFCH 

STA LASTCAL ;GOOD CAL 

CALL _PON 

CALL PRTIME 

FPRN eSl.BSlPT 

CALL LPOFF 

CALL CMPMCA ; COMPUTE MCA/WT 

CALL ' LPON 

FPRN y>2 ,V2PT 
CALL PRTHYP 

CALL LPOFF 
' MUI A.OFFH 

STA GCODBASE (GOOD BASELINE 

-•-■9RA-A-. - j 0000 CALS 
ELSE 

. FFOA 1 

FSTA RCCAL 

FSTA ABCAL 

CALL CMPMCA ; COMPUTE MCA/VT 

CALL LPON 
LXI H, STARS 

CALL TXTYP 

CALL PRTIME 

FPRN BSBAD.BSBPT 

FPRN V2.V2PT 
CALL PRTHYP 
LXI H, STARS 

CALL TXTYP 

CALL LPOFF 
M^ I A,0 

STA GOODBASE ;8AD BASELINE 
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RLE: QDCAL.ASM 

/ 

y 





;BAD CALS f 
: ISO QUANT WITH ^ASELINE 



STA CALTYPE 
CALL LPON 

CALL PRTIME ^ . " 

=?RN 533.3S3PT 
CALL LPOF? 
;CALS DO NOT CHANGE 

FFOA 1 
FSTA RCCAL 
FSTA A3CAL 

CALL CMPMCA. i COMPUTE MCA AT 

CALL LPON 
FPRN V2 ,V2PT 

CALL PRTHYP 

CALL LPOFF 
MVI A.OF=H 

STA GOODBASE i GOOD BASELINE 

ORA A i-GOOD CALS 

ELSE 

IF ( .A ,£G , 7*> STVPEO IN CALS , QUANT , JUS" GETTING BASELINE STUFF 
LOA QUALFLG 
IF <.A,IS,2EPC" 
CALL "-PON 
CALL PRT1ME 
FPRN B37.3S7PT 
CALL . LPOFF 

j CALS 00 NOT CHANGE 

FFOA 1 
FSTA RCCAL 
FSTA ABCAL 

CALL CMPMCA ; COMPUTE MCA/VT 

CALL LPON 
FPRN • V2.V2FT 
CALL PRTHYP 

CALL ■ LPOFF 

MVI A.OFFH 

STA G00D8ASE ;GOOD BASELINE 

ORA A ;GOOD CALS 

ELSE {TYPED IN CALS , QUAL SO WILL CHANGE CALS 

■ 

jSAVE CALS 

jNOW MUST USE THE MEAN VT MEASURED TO SCALE THE CAL FACTORS SO 
; MEAN VT - VTM>G'10 \ 
CALL UTSAVECL 



IF <PSW,IS,NCARRY> ;THEN GOOD CALS 
CALL LPON 
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FILE: QDCAL .ASM 



CALL PRTIME 

FPRN V7.V7PT 

CALL LPOFF 

CALL CMPMCA ; COMPUTE MCA/VT 

CALL LPQN 

FPRN V2.V2PT 
CALL PRTHYP 

CALL LPOFF 
MVI A .OF B H 

3TA 50CDBASE iGOOO 3ASELINE 

ORA A jSOOO CAL5 
ELSE 

FFOA ! 

F3TA RCCAL 

F3TA A8CAL 

CALL CMPMCA i COMPUTE MCA/VT 

CALL LPON 

CALL PRTIME 

FPRN SS8AD.3SBPT 

FPRN V2.V2PT 

call prthyp 
call . 3 -:ff 

MVI A,0 

STA G000BASE :8A0 BASELINE 

STC i BAD CAL3 

END IF 

;** 

END IF 
END IF 
END IF 

ENDIF 
RET 



CMPMCA : 



COMPUTE MEAN MCA/VT 

MUST REMEMBER THAT ALL RC VALUES MUST 8E MULT BY RCCAL 
i AND ALL AB VALUES MUST 3E MULT BY ABCAL 

NOTE: WILL NOT BE ABLE TO COMPUTE MCA/VT THE CLEAN WAY OF USING 
i THEN SUM OF THE SUM'S DERIV BECAUSE THE THE CAL FACTOR PROBLEM 

! SO WILL HAVE TO DIVIDE BY RC DELTA * RCCAL + AB DELTA * ABCAL 

UK I H.XVALUES 
SHLD ADDR1 
LDA DLTFULL 
IF ( .A.IS.NZERO) j FULL 

LHLD DLTIPTR 
SHLD OFFSET 



ft v «-r t i — 



AO 



FILE: GDCAL.ASM 



ELSE 

LXI H,0 

SHLD OFFSET 
END IF 

FFOA 0 

fsta vtsum 

fsta mcasum 

fsta xcnt 

LOOP 

LHLD ADDR1 

LDEO MCAADDR 

CALL TOHE 
IF CPSU.IS.ZEP.CO 

RET 
ENDIF 

;GET MCA VALUE 

LHLO ADDR1 

MOV E,M 

I NX H 

MOV 0,M 

SDED . RCMCA 

INX H 

MOV E,M 

INX H 

MOV D,M 

30ED ABMCA 

INX H 

MOV E.M 

INX H 

MOV 0,M 

SOED SMMCA 

INX H 

SHLD ADDR1 

FILA -RCMCA 
FLOT 

FLDB RCCAL 

FMUL 

FPSH 

FILA ABMCA 
FLOT 

FLDB A8CAL 

FMUL 

FATB 

FPOP 

FADD 

FATB 

FFOA 15 J8ECAUSE OF 9 PT DERIV 40 

5/4 BECAUSE OF DIVISION WHEN COLLECTED 

FDIV 

FLDB CCS 
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FMUL 

FSTA MCAVALUS 

i 

;GET VT VALUE 



LHLO 


OFFSET 




v i RCDhTA 


OAD 


0 


MOV 


E.M 


INX 


H 


MOV 


0 iM 


30EO 


FRC 


LHLO 


OFFSET 


LXI 


0,A80ATA 


OAD 


0 


MOV 


E.M 


IKX 


H 


MOV 


O.M 


SOEO 


FAS 


FILA 


FRC 


PLOT 




FLOB 


RCCAL 


FMUL 




FPSH 




FILA 


FAB 


FLOT 




FLOS 


A8CAL- 


FMUL 




FATS 




FPOP 




FADD 




FLOS 


CCS 


FMUL 




FSTA 


VTVALUE 



; MCA SUM = SUM OF M CAVA L U E. ''/ r VA L U E 

* 

FLOB MCAVALUE 
FLDA VTVALUE 
FOIV 

;IF < t , THEN MAKE - 1 
FPSH 
FATB 

FFOA 1 
FSUB 

FT ST AREG , LT 
IF <PSU,IS,ZERO) 

FPOP 

FFOA 1 
ELSE • 

FPOP . 
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FILE; QDCAL.ASfl 



ENDIF 



FLD8 
FADD 
FSTA 



MCASUM 



MCASIJM 



FLDB LTVALUE 

FLDA UTS'JM. 
FADD 

FSTA VTSUM 

FLD8 XCNT 

FFOA 1 
FADD 

FSTA XCNT 

* 

LHLD OFFSET 

I NX H 

inx h 

shld ofpset 

lded oltiptf. 

call tohe 
if (psw.is.zero) 

RET 
ENDI F 

•NOW OFFSET POINTS TO EXP DELTA— SKI? IT '! 

LHLD OFFSET 

I NX H 

INX H 

SHLD OFFSET 

LDED DLTIPTR 

CALL TDHE 
IF CPSW.IS.ZERO) 

RET 
ENDIF 
ENDLOOP ' 



FLDB 
FLDA 
FDIV 
FSTA 



UTSUM 
XCNT 



VTVALUE iMEAN VT 



FLDB 
FLDA 
FDIM 
FSTA 



MCASUM 
XCNT 



MCAVT 



;MEAN MCA/VT 



MEAN MCA/VT = MEAN MCA / MEAN UT 
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FLDB MCAVALUE 
FLDA VTVALUE 
FDIV 

;IF . 1 , THEN SET * i 
FPSH 
FATB 

FFOA 1 
F3UB 

—ST AREG.LT 
IF (PSU.I3.ZER0) 

FPOP 

FPQA 1 
ELSE 

FPOP 
END IF 

F5TA MCAVT 



COMPUTE THE MEAN VALUE OF •< TIME tq PEAK INSP FLOW 
AND MEAN VALUE OF *. TIME TO PEAK EXP FLOW 



LHLO 


TMPTR 


LXI 


D.DLTIMES 


CALL 


TOHE 


IF (psu,is,:erq> 


LXI 


H.O 


5HL0 


PCNUMBER3 


RET 
END IF 




FFOA 


0 


FSTA 


PKISUM 


F3TA 


PKESUM 


FSTA 


TISUM 


LXI 


H,0 


3HLD 


PCNUMBERS 


LXI 


H.DLTIMES 


SHLD 


A0DR1 



LOOP 

iMAKE SURE THAT ADDR1 + 10< TMPTR 
LHLO ADDR1 
LXI 0,10 
DAD D 
LDED TMPTR 
ORA A 
DSBC D 

IF XPSW , IS ,NCARRY> sTHEN ALL DATA NOT THERE 

RET 
END IF 

» 

i ASSUME THAT AD0R1 PT3 TO MINI 
LHLD ADDR1 
MOV E ,M 
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I NX H 

MOV 0 ,M 

SDED 7NIN1 
j GET ?KE. MAX, PKI, MIN2 VALUES 

INX H . 

MOW E,M' 

INX H 

MOW D,M 

SOED PKE ;.=-:R -AST 3REATH 

INX H 

MOV E.M 

INX H 

MOV D ,M 

SDED TMAX 

INX H 

MOV E.M 

INX H 

MOV D,M 

SOED PKI 

INX H 

SHLO AODR! :S0 NEXT TIME ADDRl PT3 TO MINI OF NEXT 

MOV E.M 

INX H 

MOV D.M 

SOED 7MIN2 

I MX H 

MOV E ,M 

INX H 

MOV D ,M 

SOED PKE ;FOR THIS BREATH 

','/. T TO PEAK INSP FLOW = i PKI -TMIN1 1 •' -MA/- T MIN1 > * lOQV; 

LHLD TMAX 

LDED TMIN1 

ORA A 

DS8C D 

SHLO TTEMP 

LHLD PKI 

LDED TMINl 

ORA A 

DS8C D 

SHLD TT2EMP 

FILA TT2EMP 
PLOT 
FAT 8 

FILA TTEMP 

PLOT 

FDIV 

FATB 

FFOri 100 
FMUL 

i'A T TO PEAK INSP FLOW 

FLDB PKI SUM 
FADD 
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FSTA PKI3UM 

-,'/. T TO PEAK EXP FLOW:* (PKE-TMAX) / (TM1N2-TMAX) « 100 '/. 

LHLD PKE 

LCED 7M.AX 
ORA 

:ssc 0 

SHLD TTEMP 

LHLD TMIN2 

LOED 7MAX 

jPA a 

033C 0 

SHlD T72EMP 

FILA 7~EMP 
=L07 
FATB 

FILA T T2SMP 
FLOT 

FAT 3 

FrOA 100 
• FMUL 

\'A T 70 PEAK EXP FLOW 

=l:b pke sum 

FADD 

=3"A PKE3UM 

i 

;TI = TMAX - TMIN1 / 3FREQ 

lHLO TMAX 

lDED TMIN1 

ORA A 

038C 0 

SHLO TTEMP 

FILA TTEMP 
PLOT 
FATB 

FILA 3FREQ 
FLOT 
FOIL* 
;AREG = TI* 

FLD8 TISUM 
FADD 

FSTA . TISUM 



LHLO PCN UMBERS 

I NX H 

SHLO PCN UMBERS 



PKISUM 
PCNUMBEP.S 



ENOLOOP 

t 
I 

FLD8 
FILA 
FLOT 
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FDIV 
F3TA 

FL08 
FILA 
PLOT 
FOIV 
F's T A 

FLOS 
p I LA 
PLOT 
FDIV 
FSTA 



PKIAVG 

PKESUM 
PCN UMBERS 



PKEAVG 

TI5UM 
PCNUM8EP3 



T I CONTROL 



jVTT I CONTROL = ONEHUNO 
-PW8- 

FLDA T I CONTROL 



\ 



I 



FDIV 
F3TA 



VTT! CONT 5 C , L 



<<? C 



<3 



NOTE: THIS COMPUTATION MUST SE DONE AFTER MCA CC J " s, JTAT!ON 
BECAUSE HYPOPV «N0 MCAVALUE ARE AT THE SFtf -DOR 
TO C0N3ERV RrtM 1 ! " 

UILL WANT TO PRINT OUT VALUE OF HYPOPNEA CVTLIMI" • = ML 
LOA QUALFL3 

IF ( .A, IS, ZERO) j QUANT 30 IN ML 
FLOA VTLIMIT 

F3TA HYPOPV jHYPOP IN ML 
ELSE ;QUAL 30 PRINT IN TERMS OF V. 

FLOS VTMX 
F=OA 10 
FDIV 

iAREG = CONTROL VALUE IN ML 

FLD8 VTLIMIT ;HYPOP VOLUME IN ML 

FDIV 

FATB 

FFOA 100 
FMUL 

FSTA HYPOPV ;HYPOP IN OF CONTROL 
ENOIF 



RET 



STARS: .ASCI 3 ' 
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8SBA0 .' .ASCI 3 ' 

*■* SAO ISQWOLUME SCALING CALIBRATION: 
TRY ED TO SCALE TO 303 ML 

RETRY WITH Y-SCALING AT LOWER VALUE ' 
333PT: .yQRO ON EH UNO 

• 
t 

? 

I 

8ADC: .A3CIS 

«** BAD QUALITATIVE CALI5RA t IjN : 
RC 5TD = 303 
hB STD = 303 
» DELTAS = 301 
MEAN RC - 303 
MEAN AB = 303 

3ADPT: 

.WORD SC3TD , ABSTD ,aCNT . FP.C . FA 8 



351 : .ASCIS ' 
BASELINE: 

QUALITATIVE I3QV0LUME CALIBRATION: 
From 301 de! tas: 

RC CAL = 303 
A8 CAL = 3D3 
Breathing. Freq = 303 Bths/min 



8S3: .ASCIS ' 

BASELINE WITH QUANTITATIVE ISOVQLUME CALIBRATION: 
8reathing Freq = 303 Bths/min 

333PT: 

.WORD FREQCONTRL 

I 

BS7? .ASCIS ' 

BASELINE WITH GAINS TYPED IN QUANTITATIVELY 
Breathing Freq = 303 Bths/min 
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8S7PT i 

.WORD FREQCONTRL 

y 

VI i .ASCIS ' 
QUALITATIVE CALIBRATION : 
"From 301 deltas: 

RC CAL = 303 
AB CAL = 303 
3reath?ng Freq = 303 3ths/min 



V7i .ASCIS ' 
3ASELINE WITH GAINS TYPED IN QUALITATIVELY i 
From 301 del tas: 

RC CAL = 303 
AS CAL = 303 
Breathing Freq ■ 303 Bths/min 



V2: .A3CI3 ' 

From 3D! breaths; 
Mean MCA/VT = 3D3 

Mean XT ?K INSP FL = 303 X 
Mean XT ?K £. - -L = 303 X 
Mean TI = 303 sec' 



HYP: .ASCIS ' 

Hypopnea volume = 303 ' 
AHYP: .ASCIS 'ml 



RHYP 2 .ASCIS 'X of basel ine 
HYPPT: .WORD HYPQPV 



V7PT: 

V1PT: .WORD XCNT .RCCAL ,ABCAL 
.WORD FREQCONTRL 

t 
« 

V2PT: .WORD XCNT 
•WORD MCAVT 
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.WORD PKIAVG.PKEAVG 
.WORD T I CONTROL 



RESCrtLSi 

"LD3 'J7MX 
FFCA 10 
FDIV 

?S7A CNBHUNO i'HIS THE THE VOLUME SCALING 70 
iCGMPUTE THE FRE3 SO CAN COMPUTE THE ','MIN FDR CONTROL 
FILA SECONDS 
PLOT 
FAT3 

FFCA 60 
FDIU 
FAT 8 

FILA ELFMIN 
FLOT 

FADO iAREG = tt MINUTES FOR COLLECTION 
FPSH 

i ■ 

LHLO MAXCNT ;H OF DELTAS 
SRAR H 

RARR L . iMAXCNT/2 = ft BREATHS 
3HL0 ITEMP 
FILA ITEMP 
FLOT 

FAT3 ;ft 3REATHS 
F°OP ;tt MINUTES 
FOIV 

FSTA FREQCONTRL 

* 

FLDB ONEHUND ; 100 V. VOLUME 
FLDA FREQCONTRL ;FREQ 
FMUL 

FSTA VMINCCNTRL ;VMIN OF. CONTROL 



vSMMSATI * MEAN VT IN CUsN 
-f=M SMH E A N y 




CCS 

MEAN [MEAN VT IN ML 



LDA QUALFLG 

IF <.A,IS,N2ER0J jQUALITATIVE 



i 
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• SO INIT VTLIMIT TO ONEH NO * VTPCLX 

FLD8 VTPCL 

FPOA 100 
FDIV 

FLD8 CNEHUND 
FMUL 

FSTA VTLIMIT ;FvR HYPOPNEA DETECTION 

• 

ELSE ; QUANTITATIVE. USE MEAN VT INSTEAD OF ON EH UNO 

FL08 VTPCL 

FFOA 100 

For; 

FLOS MEAN sMEAN VT IN ML 
FMUL 

FSTA VTLIMIT i?GR HYPOPNEA DETECTION 

FLD8 VTMX 

FFOA 1 0 
FDIV 

FL08 MEAN 
FSWP 

FDIV ; RAT 10 OF •:VTMX/tO> / MEAN VT 

FSTA MEAN 
iRCCAL = CURRENT RC GAIN >'. RCGNF> * MEAN 
jABCAL = CURRENT AB GAIN • nBGNF) * MEAN 



FLDA 


MEAN 










FLDB 


RCGNF 


iCAL 


FROM 


INITIAL 


CALCS 


FMUL 












FSTA 


RCCAL 










FLDA 


MEAN 










FLDB 


A8GNF 


?CAL 


FROM 


INITIAL 


CALC3 


FMUL 












FSTA 


A8CAL 











NOW THE CAL FACTORS ARE SCALED SO VTMX/ 10 WILL APPEAR AS :00 V. DURING 
LONG TERM 
RET 



•EXTERN QD2 

.LOC QD2 

RCMCAs . BLKB 2 

ABMCA; .BLKB 2 

SMMCAi , 8LK8 2 

ADDR1 : . 8LKB 2 

M CAVA LUES 
HYPOPVt 

. 3LKB 4 

i 
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VTWALUEt . BLKB 4 

MCAVT: .BLKB 4 

yTSUM : .BLKB 4 

MCASUM: .BLKB 4 

PKIAVGi .BLKB 4 

PKEAVG: .BLKB 4 

71 SUM I .BLK8 4 
.PELOC 

.END 
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. I DENT VSETUP 
. INSERT SP80.ASM 
.INSERT FPMAC.SRC 
. INSERT CALCOM.ASM 
. INSERT CKCOM 
. INSERT SELCOM 
.INSERT BAGCOM 
.INSERT VARCOM.CCC 
.INSERT SYCOM 
.INSERT HRDCOM 

. EXTERN AP8P,2ALL,AUOEN 

.EXTERN RTON.RTOF.RTICK 
.EXTERN TDHE 
.EXTERN DERIV5 

.EXTERN PTRUPD f 6ETP2,FIND,R0UNO 
.EXTERN SPON.SPOF 



.ENTRY SETUP 



CR = ODH 
FINIT 

; 

SETUP: 

Jsetup The clock for validation 

■ 

•SAMPLE AT 20 PTS PER SEC 

PUSH H 
; INIT REAL-TIME BUFFER POINTERS 
i MODIFIED 09-24-81 



mv I 


A, 3 


STA 


WAITFLS 


XRA 


A 


STA 


URAPF 


STA 


SECFLQ 


CMA 




STA 


8FLG 


XRA 


A 


STA 


VFRSTMX 


LXI 


H,0 


SHLD 


ELPTM 


CHLD 


ELPMIN 


SHLD 


SECONDS 


SHLD 


CPPTR 



t 



1 

! 
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LXI 


H,0 




SHLD 


RIPT 




SHLD 


ODRIPT 




SHLD 


SPTEMP 




SHLD 


PNSUM 




LXI 


H,8 


JSKIP THE 1ST 4 DATA 


SHLD , 


ROPT 




LXI 


H,TM8UFF 




mi 


M,0 




I NX 


H 




MVI 


M,0 




LXI 


H ,0 




SHLD 


TIME 




VARIABLES FOR PICKING MAX AND MINS 


XRA 


A 




STA 


PTF 


5 START BY LOOKING FOR 


STA 


SLOPEF 


J IN IT MAX NOT FOUND 


XRA 


A 




STA 


INTFLG 




CMA 






STA 


GFLG 


i INDICATES CLOCK ON 


CALL 


AUDEN 




CALL 


AP8P 




CALL 


2ALL 




LXI 


H.IRSRV 




mi 


A, 20 


(CLOCK RATE 


CALL 


RTON 


;TURN ON CLOCK' 


LDA 


BAG 




IF (.A, 


IS, ZERO) 


jUSING SPIROMETER 


CALL 


SPON 




ENDIF 






POP 


H 





RET 



CLKSTP: 



;STOP DATA COLLECTION 
iTURN CLOCK OFF 



MODIFIED 09-24-81 



CALL RTOF ;TURN OFF CLOCK 
LDA SAG 
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IF < .A, IS, ZERO) ;USING SPIROMETER 

CALL SPOF 
ENOIF 

XRA A 

STA GFLG 

CALL AUDEN 

CALL APBP 

POP P9W 

POP H 

POP D 

POP 8 

RET 

i 

IRSRV: 



INTERRUPT SERVICE ROUTINE 

PUSH 8 

PUSH 0 

PUSH H 

PUSH PSW 



CALL RTICK {SERVICE INTERRUPT 
LDA INTFLG 

IF < .A.IS.NZERO) jTHEN INTERRUPT INTERRUPT 

NOP 
ELSE 

CMA 

STA INTFLG (SET FLAG. 
END IF 

CALL CLKTIK 



CALL SAMPLE 



LDA SLOPEF 

IF ( .A IS ZERO) 

CALL* GETP2 {RETURNS C IF < 5 PTS IN BUFFER 

jTHIS JUST MAKES SURE EXTRA 
{POINTS IN BUFFER 

. JC EXIT ;NOT ENOUGH POINTS IN BUFFER 



.WO 87/04332 



PCT/US87/00217 



FILE: VSETUP.ASM 



-61- 



MVI 
STA 
ENOIF 



A,OFFH 
SLOPEF 



CALL 



DERIV5 



LDA VFRSTMX 

IF k. A, IS, ZERO) ; HAVE NOT FOUND 1ST MAX 
LHLD ROPT 
CALL FIND 

IF (PSU,IS,NCARRY) ;NOT MAX 
LHLD ROPT 
LXI 8.R3IZE 
CALL PTRUPD 
SHLD ROPT 
JMP EXIT 

END IF 

jMAX FOUND 
LDA UFRSTMX 
CMA 

STA VFRSTMX 
LDA PTF 

CMA 

STA PTF 

SHLD SPTEMP ;TEMP CRIT PT 

SDED RCTEMP i CORRESPONDING RC PT 

3BCD ABTEMP ; CORRESPONDING AB PT - 

LHLD ROPT 

LXI B.RSIZE 

CALL PTRUPD 

SHLD ROPT 

JMP EXIT 



5INIT MAX FOUND 
LHLD ROPT 
CALL FIND 

UPDATE ROPT FOR NEXT TIME 

PUSH PSW 

PUSH D 

PUSH B 

PUSH H 



LHLD ROPT 

LXI B,RSIZE 

CALL PTRUPD 

SHLD ROPT 

POP H 

POP B 



ENDIF 
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POP 
POP 



0 

PSU 



i 



PINO RETURNED C SET IF CP FOUND 
AND HI HOLDS SP PT VALUE 
AND DE HOLDS RC PT VALUE 
AND BC HOLDS A8 PT VALUE 

IF < PSU, IS, CARRY) 
SHLD SPVALUE 
SDED RCVALUE 
S8CD ABVALUE 



LDA 



PTF 



IF < .A, IS, ZERO) jTHEN MAX FOUND 



XRA 


A 


cm 




3TA 


BFLG 


IX I 


D.5PMAX 


LHLD 


CPPTR . 


DAD 


D 


LDA 


BAG 


IF (.A, IS, ZERO) 


LDED 


SPVALUE 


ELSE 




LDED 


FIXCU 


END IF 




MOV 


M,E 


I NX 


H 


MOV 


M,D 


LHLD 


CPPTR 


LXI 


D.RCMAX 


DAD 


D 


LOED 


RCVALUE 


MOV 


M,E 


1NX 


H 


MOV 


M , D 


LHLD 


CPPTR 


LXI 


D.ABMAX 


DAD 


D 



i SPVALUE HOLDS SPI=CMETER VALUE 
iUSING SPIROBAG, SO FIXCU HOLDS VALUE 



IF VALIDATING AND IF USING SPIROBAG 
THEN SPVALUE HOLDS THE SUM MAX 
LDA BAG 

IF (.A,IS,N2ER0> ;THEN USING SPIROBAG 
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L1001: 



ELSE 



LDEO 


SPVALUE 


JMP 


L1001 


END IF 




LOED 


ABVALUE 


MOV 


M,E 


I NX 


H 


MOV 


M,D 


LHLD 


MAXCNT 


ir^x 


H 


SHLO 


MAXCNT 


LHLD 


CPPTR 


1NX 


H 


I NX 


H 


SHLD 


CPPTR 


LHLD 


MAXCNT 


LXI 


D , 8RTHMX 


ORA 


A 


DSBC 


0 


IF CPSW 


, IS, ZERO) 


JMP 


CLKSTP 


END I F 




LHLLJ 


CPPTR 


LXI 


C5PMIN 


DAD 


0 


LDA 


BAG 


IF C .A, 


IS, ZERO) 


LDED 


SPVALUE 


ELSE 




LXI 


0,0 


END IF 




MOV 


M,E 


"INX 


H 


MOV 


M,D 


LHLD 


CPPTR 


LXI 


D.RCMIN 


DAD 


0 


LDED 


RCVALUE 


MOV 


M,£ 


INX 


H 



j SUM VALUE 



(THEN BUFFER FULL 



; SPVALUE HOLDS SPIROMETER VALUE 
;USING SPIRQBAG.SO MIN = 0 



MOV M , 0 

LHLD CPPTR 

LXI O.ABMIN 

OAD 0 



I 
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LHLD DATA 

MOV A ,H 
CMA 

ani OFH 

MOV H,A 

MOV A,L 
CMA 

MOV L,A 

SHLD DATA 

END IF 
RET 



i 

PRNG: 

HL HOLDS THE RING INPUT PTR 
DE HOLDS ADDR OF RING 

DAD D 

LDED OATA 

MOV M,E- 

INX H 

MOV M,D . 
RET 

.END 
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. I DENT QDSTATS 
. INSERT FPMAC.SRC 
. INSERT SP80 
FINIT 

.INSERT VTCCOM 
. INSERT CALCOM 
. INSERT OLQCOM 
.INSERT CKCOM 
.INSERT MARCOM.CCC 
.INSERT SYCOM 
.INSERT QD1.ASM 



. EXTERN TOHE ..ROUND 





. ENTRY 


PRTQDC 




.ENTRY 


CMPSTATS 


PRTQDC: 






{THIS 


!S THE ENTRY 




CALL 


CMPSTATS 




RC 






FLDB 


RCSTO 




FLDA 


ABSTD 




FDIV 






FSTA 


A8CAL 




FPRN 


QDCV, QDCPT 


* 


RET 




QDCV; 


.ASCIS 




QDC = 


303: 


» DELTAS * 




MEAN 


STD COV 


RC: 


3D3 


3D3 3D3 


A6: 

• 


303 


3D3 303 


QDCPT s 


.UORO 


ABCAL.XCNT 



3D1 



.UORD 
.UORO 



FRC,RCSTD,COVRC 
FAB, ABSTD, COVAB 



NODELTS: 

NO QDC DELTAS 



.ASCIS ' 
AVAILABLE 



CMPSTATS : 

; RCDATA AND ABDATA HOLD THE DELTAS FOR RC AND AB 

i 

; COMPUTE THE MEAN AND STD OF RC, AB AND RC+AB 
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LHLD DLTIPTR 

SHLD DLTOFF ;3AVE IT INCASE IT CHANGES IN THE 
{MIDDLE OF THIS DURING LONG TERM 



LDA DLTFULL 

IF (.A, IS, ZERO) j RING NOT FULL YET 
LXI H,0 
SHLD OFFSET 
LDED DLTOFF 
CALL TDHE. 

IF <PSU,IS,ZERO) {THEN NO DELTAS 

LXI H.NODELTS 
CALL TXTYP 

src 

RET 

END IF 
ELSE 

LHLD DLTOFF 



SHLD 
END IF 



OFFSET 



sRING IS FULL 

{POINTS TO NEXT PLACE TO PUT, SO WILL BE 
;THE OLDEST VALUE 



1ST TIME THRU THE OATA COMPUTE MEAN AND STD OF JJST SUM 

2ND TIME THUR THE DATA ONLY USE DATA UHERE SUM IS WITHIN +-QDCSTD OF MEAN 
OF SUM 



FFOA 
FSTA 
FSTA 
FSTA 



0 

SMSUM 

SMXXSUM 

XCNT 



LOOP 



LHLj) 

LXI 

DAD 

MOV 

I NX 

MOV 

SDED 



OFFSET 

D , RCDATA 
D 

E, M 
H 

D,M 
IRCV 



LHLD 

LXI 

DAD 

MOV 

INX 

MOV 

SDED 



OFFSET 

D. ABDATA 
D 

E, M 
H 

D,M 
IABV 



i 

! 

i 
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CALL SUMMERS 

LHLD OFFSET 

I NX H 

INX H 

SHLD OFFSET 

LXI 0,1000 

CALL TDHE 
IF (PSU.IS.ZERO) 

LXI H,0 

SHLD OFFSET 
END IF 



! QUIT WHEN OFFSET = DLTOFF 
LHLD OFFSET 
LDEO OLTOFF 
. CALL TDHE 
IF <PSW,IS,ZERQJ 

RET ;EXIT LOOP 
END IF 

ENDLOOP 

{COMPUTE MEAN AND STD OF SUM 





FLDB 


SMSUM 




FLDA 


XCNT 




FOIV 






FSTA 


SMMEAN 




FLDB 


SMSUM 






SMSUM 




FMUL 






FATB 






FLDA 


XCNT 




FOIV 






FLDB 


SMXXSUM 




FSUB 






FATB 






FLDA 


XCNT 




FDIV 






FSQR 






FSTA 


SMSTD 


HIGH 


LIMIT = 


SMMEAN + 


LOU 


LIMIT = 


SMMEAN - 




FLDB 


QDCSTD 




FLDA 


SMSTD 




FMUL 






FLDB 


SMMEAN 




FADD 






FSTA 


3MHIGH 




FLDB 


QDCSTD 



QDC3TD*SMSTD 
QDCSTD *SMSTD 



i 

t 



WO 87/04332 



68 



PCT/US87/00217 
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FLDA 
FMUL 
FLDB 
FSUB 
F3TA 



SMSTO 

SMMEAN 

3ML0W 



NOW COMPUTE MEAN AND STO OF RC AND AB FOR ALL PAIRS WHERE 
RC+AB < SMHIGH AND RC+A3 > SMLQW 



LDA 



DLTFULL 



IF (.A,I3,2ER0) 
LXI K,0 



3HLD 
ELSE 
LHLD 

SHLD 
END IF 



0FF3ET 



DLTOFF 



OFFSET 



iRING NOT FULL YET 



5 RING IS FULL 

{POINTS TO NEXT PLACE TO PUT, SO WILL BE 
;THE OLOEST VALUE 



FFOA 
FSTA 
FSTA 
FSTA 
FSTA 
FSTA 



0 

RCSUM 

RCXXSUM 

ABSUM 

ABXXSUM 

XCNT 



LOOP 



LHLD 

LXI 

DAD 

MOV 

INX 

MOV 

SDED 

LHLD 

LXI 

DAD 

MOV 

INX 

MOV 

SDED 

CALL 

LHLD 

INX 

INX 

SHLD 

LXI 



OFFSET 

D. RCOATA 
D 

E, M 
H 

D,M 
•IRCV 

OFFSET 

D. ABDATA 
D 

E, M 

H 

D,M 
IABV 

RASUMMERS 

OFFSET 

H 

H 

OFFSET 
0,1000 
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CALL TDHE 
IF <PSU,IS,ZERO) 
LXI H,0 
SHLO OFFSET 
END IF 



iQUIT WHEN OFFSET = DLTOFF 

LHLD OFFSET 

LOED DLTOFF 

CALL TDHE 
IF <PSW,IS,2ER0) 

RET ;E<IT LOOP 
END IF 

ENDLOOP 
i 

;IF XCNT < 2, THEN NO GOOD DATA 

FLD8 XCNT 

FFOA 2 
FSUB 

FTST AREG.LT 
IF (PSU.IS.2ER0) 

LXI H.NODSLTS 

CALL TXTYP 

STC 

RET 
END IF 

; COMPUTE MEAN OF RC AND AB 

FLDB RCSUM 

FLDA XCNT 
FDIM 

FSTA FRC 

FLDB ABSUM 

FLDA XCNT 
FDIY 

FSTA FAB 

jCOTiPUTE STD OF RC AND STD OF AB 

FLDB RCSUM 
FLDA RCSUM 
FMUL 
FATB 

FLDA XCNT 
FDIV 

FLDB RCXXSUM 

FSUB 

FATB 

FLDA XCNT 

FDIU 

F3QR 
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jMAKE 



FSTA 


RCSTD 

fx W w 1 w 


PI P>R 


ARCIIM 


PI HA 


ADC) IM 










rLUn 


A UN I 


C PS T I 1 




CI PSD 
P LVo 


A QW C f 


FSUB 




FATS 




ci rsA 


vr> it 


CPk TU 




ceno 

PbUR 




FSTA 


ABSTD 


i jai nee 
VALUcb 


ML 


FLOB 


FRC 


FLDA 


CCS 


FMUL 




FSTA 


FRC 


FLDB 


FAB 


Ct Pi A 


CCS 


PMI it 




rg IH 


PAD 


FLDB 


ABSTD 


FLDA 


CCS 


FMUL 




FSTA 


ABSTD 


FLDB 


RCSTD 


FLDA 


CCS 


FMUL 




FSTA 


RCSTD 


• 

STD/MEAN 


FLDB 


RCSTD 


FLDA 


FRC 


FDIV 




FATB 




FFOA 


100 


FMUL 




FSTA 


COVRC 


FLDB 


ABSTD 


FLDA 


FAB 


FDIV 




FATB 




FFOA 


100 


FMUL 
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FILA IRCV 
PLOT 

FSTA FRC 

FILA I ABU 
FLOT 

FSTA FAB 



jDO NOT USE IF FRC+FAB > SMHIGH OR IF FRC+FAB < SMLOW 

FLDB FRC 

FLDA FAB 
FADD 
FATB 

FLDA SMHIGH 
FSUB 

FT ST AREG.GT 
IF <P3W,IS,ZER0> 



FLDB FRC 
FLDA FAB 
FADD 
FATB 

FLDA SMLOW 
FSUB 

FTST AREG , LT 
IF <PSld,IS,ZERO> 

RET 
END IF 



FLDA FRC 

FLDB RCSUM 
FADD 

FSTA RCSUM 



FLDB FRC 

FLDA FRC 
FMUL 

FLDB RCXXSUM 
FADD 

FSTA RCXXSUM 



FSTA 

ORA 

RET 



CCVAB 
A 



RASUMMERS : 



RET 
END IF 



FLDA 



FAB 
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FLOB A8SUM 
FADD 

FSTA ABSUM 

FLOB FAB 

FLOA FAB 
FMUL 

FLDB ABXX3UM 
FAOD. 

FSTA ABXXSUM 



FLDB XCNT 
FFOA . 1 
FADD 

FSTA XCNT 
RET 



SUMMERS s 

FILA IRCV 

FLOT 

FATB 

FILA IABV 

FLOT 

FADD 

FSTA FSUM 

FLDA FSUM 
FLDB SMSUM 
FADD 

FSTA SMSUM 

i 

FLDA FSUM 
FLD8 FSUM 
FMUL 

FLDB SMXXSUM 
FADD 

FSTA SMXXSUM 

* 
I 

FLDB XCNT 
FFOA 1 
FADD 

FSTA XCNT 
RET 
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FLLEi QDCSTATS.ASM 



i 
! 

« 

.END 
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. I DENT CFIND 

. INSERT SP80.ASM 

.INSERT CALCOM 

. INSERT SELCOM 

. INSERT EFFCOM 

. INSERT DRVCOM 

. INSERT SIOCOM 

.INSERT FPMAC.SRC 

.INSERT PKCOM.ASM 

. ENTRY FIND ,GETP2 , DERI MS 

. EXTERN TDHE.PTRUPD 

FINIT 

FIND ! 

•HL = OFFSET OF CURRENT PT ON ENTRANCE 

• ON EXIT: IF CP FOUND. THEN CARRY SET AND 
} HL = SPIROMETER PT 

; DE = RIB CAGE PT 



BP 


HDUUI ICM r 1 


w ~ L. J 


rnpnee 




rnoncc 


1 VT 
LAI 


u < ortjUrr 


DAD 


Q 


MOW 


E,M 


I NX 


H 


MOV 


D,M 


SDED 


CURSP 


LHLD 


CUROFF 


LXI 


D.RC8UFF 


DAD 


D 


MOV 


E,M 


I NX 


H 


MOV 


0,M 


SDED 


CURRC 


LHLD 


CUROFF 


LXI 


D.ABBUFF 


DAD 


D 


MOV 


E,M 


INX 


H 


MOV 


D,M 


SDED 


CURAB 


LHLD 


CUROFF 


LXI 


D,TMBUFF 


DAD 


D 


MOV 


E,M 


INX 


H 


MOV 


D,M 
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SDEO 



CURTM 



CALL GTSItiN ;CHECK SIGNS OF DERIVATIVES 



LDA PTF 
IF < .A.IS.N2ER0) 
LHLO SPTEMP 



LDED CURS? 




ORA A 




DS8C D 




IF >'. PSU , 1 3 .NCARRY > ;: 


LHLD 




SHLO 


SrTSMP 


LHLO 


CURRC 


SHLO 


RCTEMP 


LHLO 


CURAB 


SHLO 


A8TEMP 


LHLD 


CURQFF 


SHLO 


OFFTEMP 


LHLO 


CURTM 


SHLO 


TMTEMP 



iLOQK FOR MIN 



SNEU TMP MIN FOUND 



rl. ADD CURRENT RCDERM PT TO RC COUNTS 
;2. ADO RCTINT TO RC COUNTS 
;3. RCTCNT <= 0 

;4. ADD CURRENT A8DERV TO AB COUNTS 
J5. AOO ABTCNT TO AB COUNTS 
[6. A6TCNT <== 0 



CALL 



UP0T1 



ORA 
RET 
ELSE 



: POSSIBLE CP SO CHECK SLOPE BETWEEN 
{CURRNT PT AND TEMP PT 
CURSP 
SPTEMP 
A 
0 

;NOCd HL HOLDS DIFFERENCE 
LDED MINCC 
ORA A 
D3BC D 

IF <PSU, IS, NCARRY) jTHEN SPTEMP IS CP 



LHLD 
LDED 
ORA 
DSBC 



LHLO 


OFFTEMP 


3HLD 


CPOFF 


LHLD 


SPTEMP 


PUSH 


H 


LHLO 


RCTEMP 


PUSH 


H 


LHLO 


• A8TEMP 


PUSH 


H 
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LHLD 

SHLD 

LHLO 

SHLD 

LHLD 

SHLD 

LHLD 

SHLD 
LHLD 
SHLD 

LHLD 
SHLD 

| 1 . STORE RC COUNTS 
;2. RC COUNTS <== RCTCMT 
,•3. RCTCNT <= 0 
;4. STORE AB COUNTS 
{5. A8 COUNTS <== ABTCNT 
;6. ABTCNT <= 0 

CALL l'PDT2 

I 

POP 3 jAB VALUE 

POP D ;RC 'VALUE 

POP H sSP VALUE 

3TC ;5ET CARRY 
PET 
ELSE 

,'NOT A LARGE ENOUGH SLOPE 

; I . ADD RCOERV TO RCTCNT 
?2. ADD ABDERV TO ABTCNT 

CALL UPDT3 
ORA A 
RET 
ENDIF 
ENDIF 

ELSE iLQOK FOR MAX 

LHLD SPTEMP 

LDED CURSP 
XCHG 
ORA A 
DSBC D 

IF CPSU.IS.NCARRY) j jTHEN CURRENT >= SPTEMP 
LHLD CURSP 
SHLD SPTEMP 
LHLD CURRC 
SHLD RCTEMP 
LHLD CURAB 
SHLD ABTEMP 
LHLD CUROFF 
SHLD OFFTEMP 



CURSP 

SPTEMP 

CURRC 

RCTEMP 

CURAB 

ABTEMP 

CUROFF 

OFFTEMP 

TMTEMP 

CPTIME sTIME OF THIS CP 
CVR'M 

TMTEMP ;NEW TEMP 
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LHLD 
SHLD 



CURTM 
TMTEHP 



CALL UP0T1 



ORA A 
RET 



ELSE ; POSSIBLE CP 
LHLD SPTEMP 
LDEO C'jRSP 
ORA A 

DS3C 0 



SNOW HL = DIFF 

lded -uric: 

ORA A 

0S3C 0 



IF < PSW , 1 3 .NCrtPP ■« -i I THEN SPTEMP IS CP 

LHLD :=--zW 

:HLO CPOFF 

LHLD l-'lf'- 

PUSH H 

uHLD RC T EMP 
PUSH 

LHLD ~S~E?' ,D 
PUSH 



SHLD CP-IME jTIME OF CP 



LHLD I'JPSP 

SHLD SPTEMP 

LHLD CURRC 

SHLD RCTEMP 

LHLD CUP.AB 

SHLD A8TEMP 

LHLD CUROFF 

SHLD 0FF7EMP 

LHLD CURTM 

SHLD TMTEMP 

CALL UPDT2 

POP 8 ;A8 VALUE 

POP 0 ;RC VALUE 

POP H :3P VALUE 



LHLD 



STC 
RET 



ELSE 



:N0 T 
UPDT3 



LARGE ENOUGH SLOPE 



CALL 

ORA 

RET 



END IF 
END IF 
END IF 



i 
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QETP2: 

SEE IF AT LEAST 5 POINTS IN RING BUFFER 

LHLD RIPT 
U.l 0,10 
ORA A 
DSBC D 

IF <PSW,I3,CARRY< ; THEN. < 5 PTS 
STC 



. ADO CURRENT ABS<PCOERM) PT TO RC COUNTS 

2. ADD RCTCNT TO RC COUNTS 

3. RCTCNT <== 0 

4. ADD CURRENT ABS(A80ERV> TO AB COUNTS 

5. ADD A8TCNT TO AS COUNTS 

6. A6TCNT <*= 0 



******* 

NOTE : DO NOT SEPEATAE INPHASE AND OUT ".OF =4A3S -Z? ~, ' 




A 



END IF 
RET 



UP0T1 ! 



LHLD 
LDED 
DAD 
SHLD 



ABSRCD 

RCIN 

D 

RCIN 



LHLD 
LDED 
DAD 
SHLD 



RCIN 

RCTICNT 

D 

RCIN 



LXI 
SHLD 



H.O 

RCTICNT 



LHLD 
LDED 
DAD 



A8SABD 
A BIN 

0 
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SHLD ABIN 



1 

LHLO ABIN. 

LDED ABTICNT 

DAD D 

SHLD ABIN 

LXI H.O 

SHLD ABTICNT 

1 

;SUM: 00 NOT SEPERATE INPHASE AND OUT OF PHASE 'VALUES 

s ADD THEM ALL 

? 

LHLD A6SSMD 

LDED SMIO 

DAD D 

SHLD SMIO 

i 



LHLO SMIO 

LDED SMTCNT 

DAD D 

SHLD SMTCNT 



LXI H,0 
SHLD . SMTCNT 



FOR PEAK FLOW STUFF: 

SMDERV = CURRENT FLOW, ABSSMD = A33<3M0ERV> 

IF ON INSP SIDE:, LOOKING FOR MAX <PTF=0) 
LDA RTF 

IF (.A, IS, ZERO) ;0N INSP SIDE 
;IF CURRENT FLOW > 0, OR IF ABSSMD = SMDERV 
LHLD ABSSMD 
LDED ' SMDERV 
CALL TDHE 

IF <PSW,IS,2ER0) {THEN CURRENT FLOW > 0 
:IF CURRENT FLOW > TEMPK 



LHLD 


TEMPK 


LDED 


SMDERV 


ORA 


A 


DSBC 


D 


IF (PSW, 1 3, CARRY) 




LHLD SMDERV 




SHLD TEMPK 




LHLD CURTM 




SHLD PKTIME 



i CURRENT FLOW > TEMPK 



END IF 
END IF 

5 1 F TTPOS > TEMPK 
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LHLO TEMPK 
LDED TTP03 
ORA A 

0S8C D 

IF <PSW, IS, CARRY) i TTPOS > TEMPK 



LHLD 
SHLD 
LHLD 
SHLD 
END IF 

LXI H,0 

SHLD 

SHLD 

LHLD 

SHLD 

SHLD 



TTPOS 
TEMPK 
TTPTM 
PKTIME 



TTPC3 
TTNEQ 
CUR7M 
TTPTM 
TTNTM 



ELSE 



SON THE EXP SIDE (LOOKING FOR MIN 



< 0, 



IF CURRENT FLOW 
LHLD 

LDED SMDERV 
CALL TDHE 
IF (PSW.IS.NZEPO) 
IF ABS< CURRENT FLOW) > TEMPK 

LHLD TEMPK 

LDED ABSSMD 

ORA A 

DS8C D 

IF (PSU, IS, CARRY) 



OR IF ABSSMD NOT= 3MDERV 
A3SSMD 



I THEN CURRENT FLOW 



;ABS< CURRENT FLOW) 



LHLD 
SHLD 
LHLD 
SHLD 



ABSSMD 
TEMPK 
CURTM 
PKTIME 



! I F TTNEG 



END IF 
END IF 
> TEMPK 

•LHLD TEMPK 
LDED TTNEG 
ORA A 

DS8C D • 

IF (PSW, IS, CARRY) j TTNEG > TEMPK 



LHLD 
SHLD 
LHLD 
SHLD 
END I F 

LXI H,0 
SHLD 
SHLD 
LHLD 



TTNEG 
TEMPK 
TTNTM 
PKTIME 



TTPOS 
TTNEG 
CURTM 
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SHLD TTPTM 
SHLD TTNTM 

« 

END IF 

RET 



UPDT2: 

STORE RC COUNTS 
5 2. RC COUNTS <== RCTCNT 
;3. RCTCNT <— 0 
;4. STORE A8 COUNTS 
;5. A8 COUNTS <== A8TCNT 
;6. A8TCNT <= 0 

LHLO RC!N 
SHLD RCINVOL 



LHLO RCTICNT 
SHLD RCIN 



LXI H,0 

SHLD RCTICNT 

LHLD AS IN 

SHLD ABINVQL 



LHLD A8TICNT 
SHLD AS IN 



LXI H,0 

shl:- abticnt 

I 

iSUMs DO NOT SEPERATE INPHASE AND OUT OF PHASE 
LHLD SMIO 
SHLD SMVOL 

! 

; 

LHLD SMTCNT 
SHLD SMIO 

i 

i 

LXI H,0 
SHLD SMTCNT 

;FOR PEAK FLOW STUFF : 

jSMDERV = CURRENT F_: ::MD = A8S<SMDERV) 
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jIF ON INSP SIDE: , LOOKING FOR MAX (PTF»Q> 
LDA PTF 

IF < .A, IS, ZERO! 1 j ON INSP 3I0E 
;TEMPK !S THE PEAK FLOW, AND PKTIME IS THE TIME OF IT 
[STORE THE TIME OF THE PEAK FLOW 
CALL PKSTORE 

LHLO TTNEG :NOW WILL 3E ON EXP SIDE 

SHLO TEMPK 

LHLO 7TN7M 

SHLD PKTIME 

L<I .1,0 

SHLD TTP03 

SHLD TTNEG 

LHLD CURTM 

SHLD TTPTM 

SHLD 7TNTM 
ELSE ;0N EXP SIDE 

{TEMPK IS THE PEAK FLOW, AND PKTIME IS THE TIME OF IT 
j STORE THE. TIME OF THE PEAK FLOW 

CALL PKSTORE 

LHLD TTPOS ;NOU WILL SE ON INSP SIDE 



SHLD TEMPK 
LHLD TTPTM 
SHLD PKTIME 
LXl .4,0 
SHLD TTPOS 
SHLD TTNEG 
LHLD CURTM 
SHLD TTPTM 
SHLD TTNTM 



UPDT3: 

;1. ADD A8S<RCDERV) TO RCTCNT 
J2. ADD ABS(ABDERU) TO ABTCNT 

LHLD A8SRCD 

LDED RCTICNT 

DAD D 

SHLD RCTICNT 



END IF 



SET 



LHLD 
LDED 
DAD 
SHLD 



ABTICNT 



ABSABD 

ABTICNT 

D 
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5 SUMS 00 NOT SEPERATE INPHASE AND OUT OF PHASE 

LHLD A8SSM0 

LDED SMTCNT 

OAO 0 

SHLD SMTCNT 



t 

iFOR PEAK FLOW STUFF: 

;SMDERV = CURRENT FLOW. wBSSNO = ABS<SMDERV> 
i 

5 IF CURRENT FLOW > 0 
LHLD 3MDERV 
LDED ABSSMD 
CALL TDHE 

if (Psw,is,:ero:' sthen current flow ) o 

;IF CURRENT FLOW > TTPC3 



LHLD TTPOS 
LDED SMDERV 
ORA A 
DSBC 0 

IF <PSW, IS, CARRY) 
LHLO SMCERV 
SHLD TT?0S 
LHLD CUR™ 
SHLD TTPTM 



ELSE 

ilF ASS'I CURRENT FLOW) > TTNEG 
LHLD TTNEG 
LDED ABSSMD 
ORA A 
DSBC D 

IF <P3W, IS, CARRY) 



LHLD 
SHLD 
LHLD 
SHLD 



ABSSMD 
TTNEG 
CURTM 
TTNTM 



ENDIF 



END IF 



RET 



PKSTORE: 



LHLD 
SHLD 
RET 



PKTIME 
PKYALUE 
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GTSIGN: 



LHLD 
MOV 
AN I 
STA 

LHLD 

MOV 
AN I 
3TA 

LHLD 
MOV 
AN I 
STA 



MOV 
LDA 



RCDERV j CURRENT RC DERIV VALUE 

A , H 

30H 

RCSIGN 

ABDERV : CURRENT AB DERIV VALUE 

A.H 

30 H 

A8SIGN 

3MDERV 
A ,H 
30.H 

SMSiGN 



8. A 

RCSIGN 
IF ':.A,EQ, .8) 
XRA A 
CMA 
STA 
ELSE 
XRA 
STA 
END IF 



iTHEN RC INPHASE W SUM 



RCINPHS 
A 

RCINPHS 



LDA 
MOV 
LDA 



SMSIGN 
8 A 

ABSIGN 



IF <.A,EQ,.8> 

XRA A 

CMA 

STA 
ELSE 

XRA 

STA 
END IF 
RET 



5 THEN AB INPHASE U/ SUM 



ABINPHS 
A 

ABINPHS 



DERIV5: 

i DERIV USING 9 POINTS 
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;*******#***#*******»*** 

I 

i DERIVATIVE USING 5 POINTS 

,-SRCOFST HOLDS OFFSET TO GET Y5 FROM RAW DATA BUFFERS 

; ASSUME: RCY 1 , RCY2 , RCY 3 . RCY 4 , ABY1 ,ABY2 .A8Y3 ,A8Y4 ALREADY SET 



LXI 


D.RC9UFF 


LHLD 


SRCOFST 


DAD 


J 


MOV 


E,M 


I NX 


H 


MOV 


D.M 


SDED 


RCY 9 


LHLD 


SRCOFST 


LXI 


D.ABBUFF 


DAD 


0 


MOV 


E.M 


I NX 


H 


MOV 


D.M 


SDED 


A8Y9 


LXI 


3.R3IZS ; UPDATE FOR NEXT TIME 


LHLD 


SRCOFST 


CALL 


PTRUPD 


SHLD 


SRCOFST 



sDERIV - \«-n«4 + C3-2>*3 * '7-3)42 + vo-4) 
■ 

LHLD 
LDED 
ORA 
DSSC 
DAD 
DAD 
PUSH 

5 

LHLD 
LOED 
ORA 
DSBC 
MOV 
MOV 
DAD 
DAD 
POP 
DAD 
PUSH 

■ 

LHLD 
LDED 
ORA 
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RCY? 
RCY1 
A 
D 

H ;<RCY?-RCY1)*2 
H ;(RCY9-RCY1 -*4 

H 

RCY 8 - 
RCY 2 
A 
D 

E.L 
D.H 

H \*2 

D ;*3 

D 

D 

H 

RCY 7 
RCY 3 
A 
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DSBC 


D 


DAD 


H 


POP 


0 


DAD 


0 


PUSH 


H 


LHLD 


RCY6 


LDED 


RC>4 


CPA 




033C 




POP 


0 


DAD 


D 


3HLD 


RCDESW 


DIVIDE 


RCDERY 


LHLD 


RC3ERW 


3RAR 


H 


PARR 




3RAR 


H 


RARR 




3HL0 


RCDEW 
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INTEGRATING IT 



DERIV ■ ■:?-l>*4 * • 3-2; *3- + (7-3)*2 * '.4-4) 



LHLO 


A3Y? 


LOED 


ABY! 


ORA 


A 


DSBC 


D 


DAD 


H 


DAD 


H 


PUSH 


H 


LHLD 


A BY 8 


LDED 


A BY 2 


ORA 


A 


DSBC 


D 


MOV 


E,L 


MOV 


D.H 


DAD 


H 


DAD 


D 


POP 


D 


DAD 


D 


PUSH 


H 


LHLD 


A BY 7 


LDED 


A8Y3 


ORA 


A 


DSBC 


0 


DAD 


H 


POP 


D 


DAD 


D 


PUSH 


H 



!<A8Y9-A8Y1)*2 
KA8Y9-ABY1M4 



;*2 
5*3 
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; • 



LHLD 


A BY 4 
HO 1 O 


LDED 


HOT 4 


OP. A 


A 


D5BLr 


n 
j 


POP 


0 


DAD 


D 




ABOER'V 


: DIVIDE 


*B0ER^ 


LHLD 


ABDER'v 


3 PAR 


H 


RARR 




SRAR 


H 


RARR 


i 


SHLD 


A60ERM 
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: RESET VALUES FOR NE>~ 

LHLD RCY2 

3HLD ROY 1 

i HLD ROY 3 

■;w L D ROY 2 

LHLD ROY 4 

3HLD ?CY3 

LHLD RCY5 

3HLD ROY 4 

. 4 L D ROY d 

3 HLD RCY5 

LHLD ROY 7 

3HLD ROY 6 

LHLD ROY 3 

SHLD R0Y7 

LHLD ROY? 

SHLD ROY 8 

• 

LHLD ABY2 

SHLD A8Y1 

LHLD A BY 3 

SHLD ABY2 

LHLD A BY 4 

SHLD A BY 3 

LHLD A BY 5 

SHLD A BY 4 

LHLD A8Y4 

SHLD A BY 5 

LHLD A8Y7 

SHLD A BY 6 

LHLD . ABY8 

SHLD A8Y7 

LHLD A BY? 

SHLD ABY8 
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LHLO RCDERV 

LDED A3DERV 

0A0 0 

SHLD SMOERV 



•..COMPUTE ABSOLUTE VALUE yr 'CCSRV ^ND ABOERV, AND SMOERV 
LHLO RCDERV 
MOU A , H 
AN I 30 H 

IF < .A.I3.NZER0"' (THEN RCDERV '.. 0 



mov 




CMA 




MOV 


L.A 


MOV 


A.H 


CMA 




MOV 


H i ri 


I NX 


H 


SHLD 


-B5RCD 


:lse 




SHLD 


ABSRCO 



END IF 



LHLO A8DER'^ 

MOV A ,H 

AN I BOH 

IF C .A.IS.NZERO' ; THEN ABOERV ■: 0 

MOV A.L 
CMA 

MOV L.A 

MOV A,H 
CMA 

MOV H,A 

INX H 

SHLD ABSA80 
ELSE 

SHLD ABSABD 
ENDIF 

LHLD SMOERV 

MOV A ,H 

ANI 80H 

IF ( .A.IS.NZERO) :THEN SMDERV < 0 

MOV A , L 
CMA 

MOV L.A 

MOV A.H 
CMA 

MOV H.A 

INX H 

SHLD ABS3MD 
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ELSE 

SHLD A8SSMD 
END IF 



RET 



GETP2: 

;3EE IF AT LEAST 10 POINTS IN RING SUFFER 

LHLD RIPT 

LXI 0,20 

ORA A 

OSBC 0 

IF <PSW. IS. CARRY) jTHEN < 10 PT3 

STC 
ELSE 



SET UP INIT VALUES IN RCY 1 , RCY2 . RCY3 . RCY4 
ABY1 ,ABY2,rtBY3.ABY4,AND SRCOFST 



LXI 


H.RC8UFF 


MOW 


E.M 


I NX 


H 


MOV 


D.M 


SOEO 


RCY1 


I NX 


H 


MOV 


E,M 


I NX 


H 


MOV* 


0,M 


3DED 


RCY2 


I NX 


H 


MOV 


E,M 


I NX 


H 


MOV 


D.M 


SDEO 


RCY3 


I NX 


H 


MOV 


E |M 


INX 


H 


MOV 


O.M 


SOEO 


RCY4 
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INX H 

MOW E ,M 

INX H 

MOW O.M 

SDED RCY5 

INX H 

MOW- E,M 

INX H 

MOW O.M 

SOED RCf6 

INX M 

MOW E,M 

INX H 

MOW O.M 

SOEO RCY7 

INX H 

MOW E,M 
INX • H 

MOW 0,M 

SOED RCYS 



LXI H.ABBUFF 

MOW E,M 

INX H 

MOW D,M 

SOEO A8Y1 

INX H 

MOW E.M 

inx h 

MOW D,M 

SOED ABY2 

INX H 

MOW E,M 

INX H 

MOW 0,M 

SOEO A8Y3 

INX H 

MOW E,M 

INX H 

■*0W D,M 

SOED A3Y4 

INX H 

MOW E,M 

iro< h 
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MOV 
SDED 



D.M 
A8Y5 



INX H 

MOV E.M 

INX H 

MOV D,M 

SDEO A8Y6 

INX H 

MOW E,M 

INX H 

MOV O.M 

SDED A BY? 

INX H 

MOV E.M 

INX H 

MOV O.M 

SDED A8Y8 

LXI H,14 ;GFSET FOR THE ?TH VALUE IN BUFFERS 

SHLD SRC0F5T 

LXI H,3 {OFFSET OF MIDDLE OF DERI 1 . 

SHLD ROPT 



LXI H.O . 

SHLD IMTGRC 

ORA A . 
END IF 
RET 



. I DENT OKFIND 

•INSERT SP80.ASM 

.INSERT CALCOM 

.INSERT FPMAC.3RC 

.INSERT EFFCOM 

.INSERT DKCOM 

.INSERT DRVCOM 

.INSERT SIOCOM 



. EXTERN CFD 
.LOC CFD 



CUROFF: 
OFFTEMP: 



.BLK8 
.8LK8 
.RELOC 



.END 
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.ENTRY DKFIND.0KGETP2 
.EXTERN TDHE 
. ENTRY DERIV5 
.EXTERN PTRUPD 
FINIT 

; 
i 

: ASSUME: THAT ALL RC VALUES, A3 VALUES. AND SUM VALUES 

WILL HAVE FIRST 4 SITS = OOCOH COMING FROM 
i A-D. 

OKFIND: 

;HL = OFFSET OF CURRENT FT ON ENTRANCE 

i ON EXIT: IF CP FOUND. THEN CARRY SET AND 

5 HL = SUM PT 

; OE = RIB CAGE PT 

: BC = ABDOMEN ? T 

; TVALUE = TIME OF CP 



PUSH 


H 


PUSH 


H 


PUSH 


H 


LXI 


D.SPBUFF 


DAD 


0 


MOV 


E.M 


I NX 


H 


MOV 


D.M 


SDED 


CURSP 


POP 


H 


LXI 


0 , RC8UFF 


DAD 


D 


MOV . 


E.M 


INX 


H 


MOV 


D,M 


SDED 


CURRC 


POP 


H 


LXI 


D,AB8UFF 


DAD 


0 


MOV 


E.M 


INX 


H 


MOV 


D.M 


SDED 


CURA8 


POP 


H 


LXI 


D.TM8UFF 


DAD 


D 


MOV 


E,M 


INX 


H 


MOV 


D.M 
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SOEO CURTM 



CALL GTS I GN ; CHECK S 1 5NS OF DERIVATIVES 



IDA PTF 

IF < .A, IS, ZERO) ; LOOK FOR MIN 
LHLD SPTEMP 
LDEO CURSP 
ORA A 
DS3C 0 

IF <PSU.IS,NCAPRY> sNEU TMP MIN FOUND 

LHLD CURSP 

SHLD SPTEMP 

LHLD CUPRC 

SHLD °CTEMP 

LHLD CURA8 

SHLD ACTErtP 

LHLD CURTM 

SHLD SMTIME 



i 1 . ADD CURRENT RCDERV *! TO RC COUNTS 
:2. ADD RCTCNT TO RC COUNTS 
S3. RCTCNT <== 0 

;4. ADD CURRENT ABDERU TO AB COUNTS 
i5. ADD A8TCNT TO AB COUNTS 
;6. ABTCNT <~ 0 

CALL UPDT1 



ORA A 
RET 

ELS E. .-P0SSI8LE CP SO CHECK SLOPE 8ETLJEEN 
jCURRNT PT AND TEMP PT 

LHLD CURSP 

LDED SPTEMP 

ORA A 

DS8C D 

;NOW HL HOLDS DIFFERENCE 

LDED MINCC 

ORA A 

DSBC D 

IF (PSU,IS,NCARRY) 5 THEN SMTEMP IS CP 

LHLD SPTEMP 

PUSH H 

LHLD RCTEMP 
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Dl ICU 


u 

n 


unLU 


& DT CM D 

ho 1 tnr 


PI IQW 


u 
n 


i ui n 
LnLU 


CMT IMC 

an 1 1 PI c 


SHLD 


TVALUE 


LHLO 


CURSP 


SHLD 


3PTEMP 


LHLD 


CURRC 


SHLD 


RCTEMP 


LHLD 


CjPAB 


SHLD 


mBTEMP 


LHLD 


CURTM 


SHLD 


5MTIME 



1 . st: = = ?C counts 

2. RC COUNTS <== RCTCNT 

3. RCTCNT <== 0 

4. STORE AB COUNTS 

5. AB COUNTS <= ABTCNT 

6. A8TCNT <== 0 

CALL UPDT2 



POP 3 ;AB VALUE 

POP D i RC VALUE 

SET ~-E - S3 »0PR!ATE BIT ON RC VALUE TO INDICATE 
A MINIMUM 8RTH 

'MOV A , D 

AN I OFH jFIRST 4 BITS SHOUL 

MOV D,A 

POP H ;SM VALUE . 

STC iSET CARRY 
RET 
ELSE 

jNOT A LARGE ENOUGH SLOPE 

1 . ADD RCDERV TO RCTCNT 

2. ADD ABDERV TO ABTCNT 

CALL UPDT3 

ORA A 
RET 
END IF 
END IF 

ELSE ;LOOK FOR MAX 

LHLD SPTEMP 

LOED CURSP 
XCHG 
ORA A 
DSBC D 
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IF <PSU,IS,NCARRY) 



LHLD 


CURSP 


SHLD 


3PTEMP 


LHLD 


CURRC 


SHLD 


RCTEMP 


LHLD 


CURAB 


SHLD 


«8TEMP 


LHLD 


CURTM 


SHLD 


SMTIME 


CALL 


'JPDT1 



[THEN CURRENT >= SMTEMP 



CRA A 
?ET 

ELSE ;POSSI 
LHLD 
lDED 
ORA A 
OSSC 
{MOW HL = DIFF 
LDED 
ORA A 
OSSC 

IF ':psw,is 

LHLD 

PUSH 

LHLD 

PUSH 

LHLD 

PUSH 

LHLD 

SHLD 

LHLD 

SHLD 

LHLD 

SHLD 

LHLD 

SHLD 

LHLD . 

SHLD 



BLE CP 
, SPTEMP 
CURSP 

D 

*!NCC 



.N CARRY) 

SPTEMP 
u 

RCTEMP 
H 

ASTEMP 
H 

SMTIME 

TVALUE 

CURSP 

SPTEMP 

CURRC 

RCTEMP 

CURAB 

ABTEMP 

CURTM 

SMTIME 



iTHEN SMTEMP IS C? 



CALL 

POP 

POP 

SET THE APPROPRIATE 
A MAXIMUM 8RTH 

MOV 
ORI 
MOV 



UPDT2 

B 

0 



;AB VALUE 
jRC VALUE 



8IT ON RC VALUE TO INDICATE 
A ,D 

10H ; FIRST 4 8 ITS SHOULD 
D,A 



= 00018 



POP 



H 



iSM VALUE 
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STC 
RET 



ELSE 



CALL 
ORA 
RET 



jNOT LARGE ENOUGH SLOPE 

UPDT3 

A 



ENOIF 
ENDIF 
ENOIF 



UPDT1 : 

5l. ADO CURRENT A8S<RCDERV> PT TO RC COUNTS 
j2. ADD RCTCNT TO RC COUNTS 
;3. RCTCNT <== 0 

;4. ADD CURRENT ABS(ABDERV) TO AS COUNTS 
i5. ADD ABTCNT TO AB COUNTS 
-,6. ABTCNT ■■== 0 



LHLD A8SRCD 
LDA RCINPHS 
IF ( .A.IS,NZERO 
LDED RCIN 
DAD D 
SHLD RCIN 



ELSE 



;PC OUT OF PHASE WITH SUM 



LDED RCOUT 
DAD D 
SHLD RCOUT 



ENDIF 



LHLD 
LDED 
DAD 
SHLD 



RCIN 
RCTICNT 

D. 

RCIN 



LHLD 
LDED 
DAD 
SHLD 



RCOUT 

RCTOCNT 

D 

RCOUT 



LXI 

SHLD 

SHLD 



H,0 

RCTICNT 
RCTOCNT 



LHLD ABSABD 
LDA A8INPHS 
IF ( .A.IS.NZERO) 



;THEN AB INPHASE WITH SUM 
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LDEO MB IN 
DAD 0 
SHLD AS IN 

ELSE ;A8 OUT OF PHASE WITH SUM 

LDEO ABOUT 
DAD 0 
SHLD ABOUT 

END IF 

LHLD A3!N 

LOED ABTICNT 

DAD D 

SHLD ABIN 

LHLD ABOUT 
LDED AB T QCNT 

^ah r, 
SHLD ABOUT 

U<I H , ] 

shld ' -=":c?r 
SHLD -bt:cnt 

SUM: DO NOT SEPERATE INPHASS mND OUT OF PHASE U-..E3 
ADD THEM ALL 

LHLD A8SSMD 

LDED 3MI0 

OrtD D 

SHLD SMIO 



LHLD SMIO 

LDED SMTCNT 

DAD 0 

SHLD SMTCNT 



UK I H.O 
SHLD SMTCNT 
RET 



UPDT2: 



1 . STORE RC COUNTS 

2. RC COUNTS <== RCTCNT 

3. RCTCNT <== 0 

4. STORE AB COUNTS 

5. AB COUNTS <= ABTCNT 

6. ABTCNT <== 0 
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i til f\ 

LHLD 


RCIN 


SHLO 


RCINVOL 


1 Ut Pi 

LULU 


orn\ it 


oriLD 


RCOUTVuL 


LHLD 


RCTICNT 


SHLO 


RCIN 


LHLD 


RCTQCNT 


SHLO 


RCQlfl 


L<1 


H.O 


SHLO 


ft f^T f iT 


SHLO 


RCTOCNT 


1 Ml A 

LnLO 


AO f kl 

ABIN 


SHLO 


ABINUOL 


LHLD 


ABOUT 


oHLO 


ABQUTvQL 


LHLO 


A8TICNT 


SHLO 


A 8 IN 


LHLD 


A8T0CNT 


SHLO 


A80UT 


LXI 


H.O 


SHLO 


A8TICNT 


SHLO 


A8T0CNT 



;SUM: DO NOT SEPSRATE INPHASE AND OUT OF PHASE 
3MI0 
SMUOL 



SMTCNT 
SMIO 



H,0 

SMTCNT 



UPDT3: 

;1. ADO ABS<RCOERV> TO RCTCNT 
:2. ADD ABS(ABOERV) TO ABTCNT 



LHLD 
SHLO 

LHLD 
SHLO 

LXI 
SHLO 
RET 
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LDA RCINPHS 
IF (.A,IS,NZER0> 

LHLD ABSRCD 

LOEO RCTICNT 

DAD 0 

3HLD RCTICNT 
ELSE 

LHLD ABSRCD 

LOED RCTICNT 

DAD D 

5HLD RCTOCNT 
END IF 

10A A8INPHS 
IF < .A.I3.N2ERO) 

LHLD A8SA8D 

LDED ABTICNT 

DAD D 

3HLD ABTICNT 
ELSE 

LHLD A8SABD 

LDED ABTOCNT 

DAD D 

SHLD ABTOCNT 
END IF 



j THEN INPHASE 



DO NOT SEPERATE INPHASE AND OUT OF PHASE 

LHLD AB3SM0 

LDED 3MTCNT 

OAD 0 

SHLD SMTCNT 

RET 



GTSIGN: 



LHLD RCDERV j CURRENT RC DERIV VALUE 

MOV A , H 

AN I 80 H 

STA RCSIGN 



LHLD 
MOV 
AN I 
STA 



A8DERV 
A , H 
30H 

A8SIGN 



i CURRENT AB DERIV VALUE 



LHLD 



3MDEP.V 
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MOV 
AN I 
STA 



A,H 
80H 

SMSIGN 



MOV 
LDA 



8,A 

RCSIGN 



IF < .A.EQ, .8) 

XRA A 

CMA 

STA 
ELSE 

XRA 

STA 
END IF 



RCINPHS 
A 

RCINPHS 



LDA 
MOW 
LDA 



SMSIGN 
8, A . 
ASSIGN 



IF <.A,EQ..8.1 

XRA A 

CMA 

STA 
ELSE 

XRA 

STA 
END IF 
RET 



A8INPHS 
A 

A8INPK3 



THEN RC INPHASE W/ SUM 



;THEN AS INPHASE W SUM 



DERIV5: 

; **#*#»**«*#*###******«.♦ 
;DERIL> USING 9 POINTS 

• 

:OERIVATIUE USING 5 POINTS 

* 

i 

iSRCOFST HOLDS OFFSET TO GET Y5 FROM RAW DATA BUFFERS 

J 

ASSUME: RCYt,RCY2,RCY3.RCY4 ) A8Y1 ,A8Y2,ABY3,A8Y4 ALREADY 



LXI 

LHLD 

DAD 

MOV 
I NX 



D, RC8UFF 
SRCOFST 
D 

E, M ' 
H 
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MOV 


D,M 


SDED 

• 


RCY9 


LHLD 


3RC0F3T 


LXI 


D.ABBUFF 


DAD 


D 


MOV 


E.M 


I NX 


H 


MOV 


D.M 


SDED 

« 


ABY9 


LXI 


B.PSIZE ; UPDATE FOR NEXT TIME 


LHLD 


SRC0F3T 


CALL 


PTRUPD 


3HLD 


SPCOFST 


jDERIU = 


1>«4 + C8-2>*3 + ( 7-3) «2 + (6 


; 

LHLD 


RCY? 


LDED 


RCY1 


ORA 


A 


DSBC 


D 


DAD 


H ;* ~C V ?-SCY1>«2 


DAD 


H ; C PCY^-RCYi ') *4 


PUSH 

• 


H 


LHLD 


RCY 3 


LDED 


PCY2 


ORA 


A 


DSBC 


D 


MOV 


E>L 


MOV 


D,H 


DAD 


H ;*2 


DAD 


D ;*3 


POP 


D 


DAD 


D 


PUSH 


n 


LHLD 


RCY 7 


LDED 


RCY 3 


ORA 


A 


DSBC 


0 


DAD 


H 


POP 


D 


DAD 


D 


PUSH 


H 


LHLD 


RCY6 


LDED 


RCY 4 


ORA 


A 


DSBC 


D 


POP 


D 


DAD 


D 


SHLD 


RCOERV 



VT V/ Of /VT^* 
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;**»♦*: DIVI0E RCDERV 8Y 4 TO AVOID OVERFLOW WHEN INTEGRATING IT 



<3-2>-*3 * •: 7-3X2 + <4-4> 



DIVI0E 


RCDERV 


LHLD 


RCDERV 


SRAR 


H 


RARR 


L 


SRAR 


H 


RARR 


L 


SHLD 


RCDERV 


= >: 9-1) -*4 + C 


LHLO 


A3Y9 


LOEO 


ABY1 


ORA 


A 


DS8C 


D 


DAD 


H 


DAD 


H 


PUSH 


H 


LHLD 


A3Y3 


LDED 


ABY2 


ORA 


A 


DSBC 


0 


MO».; 


E,L 


MOV 


0,H 


DAD 


H 


nan 


u 


POP 


D 


DAD 


D 


PUSH 


H 


LHLD 


A BY 7 


LDED 


A8Y3 


ORA 


A 


DSBC 


D 


DAD 


H 


POP 


D' 


DAD. 


D 


PUSH 


H 


LHLD 


ABY6 


LOED 


A8Y4 


ORA 


A 


DSBC. 


D 


POP 


D 


DAD 


D 


SHLD 


ABDERV 


DIVIDE ABDERV 


LHLD 


ABDERV 


SRAR 


H 


RARR 


L 



;>'ri8Y?-ABYl>*2 
; 1 A BY 9 -a BY! »*4 



;*2 
i*3 



;****»: DIVIDE ABDERV BY 4 TO AVOID OVERFLOW WHEN INTEGRATING IT 
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SRAR H 
RARR L 
SHLD ABDERV 



; RESET VALUES FOR NEXT TIME 
LHLD RCY2 
SHLD RCY1 
LHLD RCY3 
SHLD RCY2 
LHLD RCY4 
SHLD RCY3 
LHLD RCY5 
SHLD RCY4 
LHLD RCY6 
SHLO RCY5 
LHLD RCY7 
SHLO RCY4 

LHLD RCYS 

SHLO RCY? 

LHLD RCY 9 

SHLD RCY 8 

i 

LHLD a8Y2 

SHLD ABYl 

LHLD A8Y3 

SHLD ABY2 

LHLD A BY 4 

SHLD A BY 3 

LHLD A8Y5 

SHLO A BY 4 

LHLD A BY 4 

SHLD A8Y5 

LHLD ABY7 

SHLD ABY6 

LHLD A8Y8 

SHLD ABY7 

LHLD ABY? 

SHLD ABY8 



LHLD RCDERV 

LDED ABDERV 

DAD 0 

SHLD SMDERV 



!C0MPUT LHLD S0LUTE " JALUE ° F RCDERU AN ° A8DERVl AND SMDEW 
MOV A , H 
AN I 80 H 

IF ( .A.IS.NZERO) [THEN RCDERV < 0 
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MOV 


A,L 


CMA 




MOV 


L,A 


MOV 


A.H 


CMA 




MOV ' 


H.A 


INX 


H 


5HLD 


A8SRCD 


ELSE' 




SHLO 


A8SRCD 


ENOIF 





MOV A.L 
CMA 

MOV L, A 

MOV A ,H 
CMA 

MOV H .A 

INX H 

SHLO A8SA8D 
ELSE 

SHLO A8SA8D 
ENOIF 

LHLD SMOERV 

MOV A.H 

AN I 30 H 

IF ( .A.IS.NZERO) sTHEN SMOERV < 0 

MOV A, L 
CMA 

MOV L,A 

MOV A ,H 
CMA 

MOV H,A 

INX H 

SHLO ABSSMO 
ELSE 

SHLO ABSSMO - 




A8DERV 
A ,H 
80H 



A.IS.NZERO) 



:7HEN A80ERV < 0 



ENOIF 



RET 
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DKGETP2: 

iSEE IF AT LEAST 10 POINTS IN RING BUFFER 

LHLD RIPT 
LXI 0,20 
ORA A 
OSSC D 

I? -PSW, IS, CARRY » [THEN < 10 PT3 

3TC 
ELSE 

■ 
l 

! 

;SET UP INIT VALUES IN RCY 1 , RCY2 » RCY3 , RCY4 
J A8Y1 ( ABY2,ABY3,A8Y4,AND 3RC0F5T 



LXI 


H.RC8UFF 


MOV 


E,M 


una 


H 


MOV 


D,M 


SOED 


RCY1 


I NX 


H 


MOV 


E.M 


I NX 


H 


MOV 


D.M 


SOED 


RCY 2 


I NX 


H 


MOV 


E.M 


I NX 


H 


MOV 


D.M 


SOED 


RCY3 


I NX 


H 


MOV 


E,M 


I NX 


H 


MOV 


D,M 


SDED 


RCY 4 


I NX 


H 


MOV 


E,M 


I NX 


H 


MOV 


D.M 


SOED 


RCY 5 


I NX 


H 


MOV 


E,M 


INK 


H 


MOV 


D.M 


SDEO 


RCY 6 
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INX 


H 


MOV 


E,M 


I MX 


H 


mqv 


D.M 




RCY7 


*liA 


u 

n 


Mfil J 
■ . (Jv 


P M 


INX 


H 


MOW 


D.M 






LA * 


~ i n O 0 U r ^ 


MQV 


E,M 


INX 


H 


ri Uv 


r> m 

U il i 




no I • 




u 

n 


MOM 


E,M 


INX 


H 


« . uv 


L> M 


SDED 




INX 


H 


MOV 


E,M 


INX 


H 


MRU 


n M 


oner* 


Q v ^ 

MO 1 %5 


1 n A 


u 
n 


MOV 


E,M 


INX 


H 








HOlt 




U 

n 


MOW 


E,M 


INX 


H 


1 1Uv 


n m 


OUCU 


MOT J 




u 

n 


MOW 


E,M 


INX 


H 


MOW 


D,M 


SDED 


A BY 6 


INX 


H 


MOW 


E|M 


INX 


H 


MOW 


D.M 



WO 87/04332 PCT/US87/00217 

t 

-107- 

FiLE: CFIND .ASM 



SDED ABY7 

I NX H 

MQY E.M 

INX H 

MOV D,M 

SDED A8Y8 

LXI H,:4 : =0* ~-E *TH '/AuJE !N SUFFERS 

S-O 3RC0F3" 

L<! H,3 ; IFFEE - *. r m ::*-.= I" - 

SHLD ?.0P T 

-XI H.O 

SH-J INTCFC 

ORA A 
END IF 
RET 
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10 Br. 

. i dent prntval 
. insert fpmac.src 
. insert sp30.asm 
. insert calcqm 
.insert lp2cqm 
. insert sycgm 

. extern prtime 

.extern lpon.lpoff 

?IN IT 

.£<"=SN ROUND 
. ENTRY PRN T '/AL 

THIS ROUTINE PRINTS ILT 'hi :ELTA v'ALUESi 

SUM SP 5UM/SP 



THERE ARE v2-«ttAXCN' 



PRNTVAL : 

i 

.•ASSUME VALUE OF FCN~ Ij IN -REG ON ENTRY TO THI • : jUTINE 
! 

F3TA FCNT 
;ONLT MAKE HrtRD copy 
SGI j 

LHLO MAXCNT 

OCX H 

SLAR L 

RALR H 

SHLD CNT 

LIQI 4 

LXI H.XUALUES 

SHLD X2PTR 

Lil H.O 

SHLD CPPTR {OFFSET FOR SP.RC.AB 
XRA A 

STA PTF ;THI3 IS ZERO WHEN 

jODD ft TIME THRU LOOP • 

XRA A 
CMA 

3TA HRDFIG 

i 

MVI A. 10 
STA LCNT 
i ONLY MAKE HARD COPY 
CALL LPON 
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j OUTPUT SEVERAL LINE FEEDS 
LXI H.CRLFS 
CALL TXTYP 
CALL PRTIME 
FPRN HEAD.HOPT 
CALL LPOFF 

L00P1 : 

!GET THE OELTA FOR SPIROMETER 

LHLO CPPTR 

LA I O.SPMAX 

DAD 0 

MOV E.M 

I NX H ' 

MOV j.M 

50ED 3 RVALUE 

LHLO CPPTR 
LDA PTF 

IF < .A,IS,NZSRO.> jTHEN EVEN * THRU LOOP 
INX H 
I NX H 

END IF 

LXI D.SPMIN 
DAD D 
MOV E ,M 
INX H 
MOM 0 ,M 

|DE = HIN SP <=T 
LHLD SPWALUE jHL = MAX SP PT 
ORA A 
0S8C 0 

3HLD SPOELT ;SP OELTA 

FILA SPOELT 

FLOT 

FL08 CCS 
FMUL 

FSTA SPOELT ;NOW I NT E RMS OF ML 
jNOU GET SUM OELTA VALUE 



LHLO CPPTR 

LXI 0 , ABMAX 

OAO 0 

MOV E.M 

INX H 

MOV 0 ,M 

SDED A8VALUE 

LHLO CPPTR 
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LDA PTF 

IF <.A,IS,NZERQ) ;THENEVEN » THRU LOOP 
I NX H 
INX - H 

END IF 

LXI 0.A8MIN 
OAO 0 
MOM g,M 
INX H 
MOV 0,M 

:0E = MIN SUM PT • 
LHLO A8VALUE ;HL = MAX SUM PT 
ORA A 
DS8C D 

'SHLO ABDELT :SUM DELTA VALUE 

FILA A80ELT 

FLOT 

FLD3 CCS 
FMUL 

FSTA ABDELT jNOU IN ML 



;NOW GET XVALUE (SUM/S») 
LDQR X2PTP. 
FOWL 

FSTA XVAL 
3DQR X2PTP 

i 

i ONLY MAKE HARD COPY 
CALL LPON 
LDA PTF 

IF •: .A.IS.2EP0) ; THEN INSP 

LXI H.I LAB 
ELSE 

LXI H .SLAB 
END IF 

CALL TXTY D 

FPRN OELT ,DELT?T 

CALL LPOFF 



LDA PTF : 

IF < .A,IS,NZERO> iTHEN EVEN tt THRU LOOP 

LHLD CPPTR 

INX H 

INX H 

3HLD CPPTR 
END IF 

LDA PTF 

CMA 

STA PTF 
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»»t - 

OSKZ CNT,L00P1 

! 

LPll: 

j CALCULATE. AND PRINT OUT THE MEAN , STO, 3T ERROR, AND V. ERROR 



FLOB 




FLDA 


FCNT 

r Will 


■ W * f 


r 3 i n 


M PAN 


FLOP 

I U U U 


a o un 


C LDA 




FMUL 




FATB 




FLDA 


FCNT 


FOIV 




FLD3 


XX SUM 


> Jww 




FAT3 




FLDA 


FPMT 
r uin i 


FDIU 




FSQR 




F3TA 


3TDY 


FLDA 


• Wi 4 i 


' j u r» 






CT IM J 






FSTA 


0 i Z. rt ?\ 


FLDB 


MEAN 


c F0A 


1 

4, 


F3UB 




FAT 8 




FFOA 


100 


FMUL 




F3TA 


PCERR 


FLDA 


PCERR 


FA8S 




FATB 




FFOA 


1 


FSUB 




FSTA 


FTEMP 


FTST 


FTEMP.LT 


IF <PSU.IS,2ER0) 



FFOA 0 
FSTA PCERR 
ENDIF 

i 

002: 
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(ONLY MAKE HARD COPY 
CALL LPON 
FPRN STATS, STPT 

iOUTPUT SEVERAL LINE FEEDS 
LXI H.CRLFS 



CALL 
CALL 



RET 



TXTYP 
LPOFF 



STATS i .ASCI 3 ' 

Mean SUM/SP = 3D3 

Stdw SUM/SP = 303 

St Err SUM/SP = 303 



% ERROR 



3D1 



STPT : .WORD MEAN , S7DU . 5TER3 , PCERR 

!LAB: .ASCI 3 '(I) 
ELAB : .ASCI 3 ''.E) 



OELTs .ASCIS *3D3 



jD3 



3D3 



DELTPT : 



.WORD ABOELT , SPDELT „<VAL 



HDPT: 



HEADs .ASCIS ' 

Val idat i on : 
SUM SP 



SUM/SP 



CRLFS: .ASCIS 



. EXTERN PRT 
.LOC PRT 
FTEMP: . BLKB 4 



SPDELT: .BLKB 

RCDELTi .BLKB 

Y2PTRi .BLKB 2 

X2PTR: .BLKB 2 

ABDELT: .BLKB 



4 
4 
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.RELOC 
.END 



1 1 3 



LCNT: 


. BLKB 


1 


FCNTi 


.BLKB 


4 


j 

MEAN: 


.BLKB 


4 


STOVs 


.BLKB 


4 


3TERR : 


.BLKB 


4 


PCERR: 


.BLKB 


4 


LQOPCNT: 


i 
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jIF VALIDATING AND IF USING SPIROBAG 
I THEN SPVALUE HOLDS THE SUM MIN 
LDA BAG 

IF < .A,IS,NZERO) iTHEN USING SPIROBAG 

LDED SPVALUE jSUM VALUE 

JMP L2001 

END IF 

LDED ABVALUE 

L2001: 

MOV M,E 
I NX H 
MOV M,0 

END IF 
LDA PTF 
CMA 

STA PTF 
ENDIF 

SKIP5i 
EXIT: 

* 

CALL CSTS 
IF <PSW,IS,NZERO) 
CALL CINP 

IF <.A,EQ,CR> i IF CARRIAGE RETURN 

JMP CLKSTP 
ENDIF 
ENDIF 

XRA A 

STA INTFLG 

« 

» 

;**»********•*»»»»»*» 
POP PSU 
POP H 
POP 0 
POP B 

RET 



CLKTIK: 



LHLD ELPTM 
INX H 
SHLD ELPTM 
LDED SFREQ 
CALL TDHE 

IF (PSU,IS,2ER0) ;THEN 1 SECON 
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LX I H,0 
SHLO ELPTM 
LHLD SECONOS 
INX H 
SHLD SECONDS 
XRA A 
CMA 

STA SECFLG 
LXI 0,40 
CALL TDHE 

IF <PSW,IS,ZERO> ;THEN 1 MINUTE 

U<I H,0 

SHLO SECONOS 
LHLD ELPMIN 
INX H 

SHLD ELPMIN 
ENDIF 
ENDIF 
RET 



SAMPLE : 

; 

LDA BAG 

IF (.A, IS, ZERO) {THEN USING SPIROMETER 

LHLD SMAD ;SUM VALUE 

NOV A , H 
CMA 

ANI OFH 

MOV H,A 

MOV A, L 
CMA 

NOV L,A 

SHLD DATA 

LHLD RIPT 

LXI D.ABBUFF 

CALL PRNG 

i 
■ 
i 

LHLD SPAD ; SPIROMETER VALUE 

MOV A ,H 

CMA 

ANI OFH 

NOV H,A 

NOV A,L 
CMA 

• MOV L,A 

SHLO DATA 

LDA PNFLAG 

IF (.A,IS,N2ER0) jTHEN INPUTTING PNEUMOTACH-SO INTEGRATE 



i 
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;NOTE:*»»* ASSUME 0 FLOW IS AT 7FFH 
LHLO DATA 
LXI D.7FFH 
ORA A 

DS8C 0 
DIVIDE BY 3 (SIGNED SHIFT RIGHT) 
CHANGE TO SHIFT ONLY 2 TIMES 
28YTES PER INSTRUCTION 

NOP 

NOP 

• 

NOP 
NOP 

', SRAR H 
; RARR L 

SRAR H 

RARR L 

SRAR H 

RARR L 

LDED PNSUM 

DAD D 

SHLD PNSUM 
SADO OFFST SO ALL * > 0 

LXI D.800H 

DAD D 

SHLD DATA 
END I F 

CALL CKFLI.P 
LHLO RIPT 
LXI D.SPBUFF 
CALL PRN6 

j 

ELSE ) USING FIXED VOLUME 

i STORE SUM DATA IN SP8UFF SO WILL 8E USED TO PICK OFF BREATHS 

» 

5 SUM VALUE 



LHLD 


SMAO 


MOV 


A,H 


CMA 




AN I 


OFH 


MOV 


H,A 


MOV 


A,L 


CMA 




MOV 


L,A 


SHLD 


DATA 


LHLD 


RIPT 


LXI 


D.SPBUFF 


CALL 


PRNG 



i'FILL ABBUFF WITH ZEROES 
LXI H,0 
SHLD DATA 
LHLD RIPT 
LXI D , ABBUFF 



i 
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CALL PRN6 



END IF 



{FILL RCBUFF WITH ZEROES 

LXI H,0 

SHLO DATA 

LHLD RIPT 

LXI 0, RCBUFF 

CALL PRNG 

i 

i 



INCREMENT RIPT FOR NEXT TIME 

LHLD TIME 

I NX H 

SHLO TIME 

SHLO DATA 

LHLD RIPT 

LXI DjTMBUFF 

CALL PRNG 

LHLD RIPT 

SHLO ODRIPT 

LXI B.RSIZE 

LHLD RIPT 

CALL PTRUPD 

SHLD RIPT 

LDA URAPF 

IF < .A, IS, ZERO) ;SEE IF RIPT RESET TO 0 

LXI D,0 

LHLD RIPT 

CALL TDHE 

IF <PSU,IS,ZERQ) 
MVI A,OFFH 
STA URAPF 

ENDIF 
END IF 
RET 



CKFLIP; 

ICHECK IF NEED TO FLIP 
LDA FLIPg 
IF ( .A,IS,NZERO) 
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5 in practice, the Respigraph ™ scaling amplifiers for the 
rib cage and abdomen signals are initially set at unity gain. \ 
10 minute baseline is then taken. The program automatically 
computes the delta values for the rib cage and abdomen signals, 
discards the wild points, computes the mean values for the rib 
15 cage and abdomen from the totals of the remaining delta values, 

computes the standard deviations for those means, and then 
2q computes the proportionality factor Z from Equation M. The 
scaling amplifier for the rib cage is then adjusted to the 
proportionality factor Z. The real time output from the 
25 Respigraph ™ representing the sum of the rib cage and abdomen 
signals is then proportional to tidal volume as explained more 
3Q fully above. The program also recomputes the proportionality 
factor Z during each subsequent 5 minute interval, if the 
recomputed value differs from 1.0 by more than an acceptable 
35 amount, the calibration routine can be rerun. The program also 
computes the scaling factor M from Equation N if instructed to 
do so upon inputting an actual tidal volume value as derived, 
e.g., by spirometry. 

While there has been shown and described herein a preferred 
45 embodiment of the present invention and certain suggested 

modifications thereto, it will be apparent that further changes 
_^ and modifications may be made without departing from the spirit 
and scope of the invention. Accordingly, the foregoing 
description should be construed as illustrative and not in the 
55 limiting sense, the scope of the invention being defined by the 
following claims. 
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Claim 1 

10 1 ' In a method for non-invasively measuring a subject's 
respiration volume of the type including providing a signal 
responsive to a rib cage dimension indicative of rib cage 

15 contribution to respiration volume, providing a signal 

responsive to an abdominal dimension indicative of abdominal 

2Q contribution to respiration volume, multiplying at least one of 
said rib cage and abdominal signals by predetermined weighting 
factors reflecting the relative contributions of said rib cage 

25 and abdomen to respiration volume, and summing said weighted 
signals for providing a signal proportional to respiration 

3q volume, the improvement comprising non-invasively determining 
said weighting factors by 

(a) totaling the delta values for said rib cage signal 

35 over a baseline period of substantially steady state 
breathing; 

40 (b) tota ling the delta values for said abdomen signal 

over a baseline period of substantially steady state 
breathing; 

45 (c) dividing the average variability of the mean of 

the total of said delta values for one of said rib cage or 

abdomen 

50 

signals by the average variability of the mean of the total 
of said delta values for the other of said rib cage or abdomen 
55 signals; and 

(d) and multiplying said other signal by a weighting 
factor equaling the quotient derived from step (c). 
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